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I.  INTRODUCTION 


A,  BACKGROUND 

In  the  recent  years,  the  United  States  Military  has  taken  an  increasing  interest  in 
developing  the  ability  to  offload  vehicles  and  machinery  from  large  roll-on  and  roll-off 
(RO/RO)  vessels  in  open  seas.  Several  classes  of  vessels  provided  for  this  purpose, 
including  are  LMSRs,  CAPE-Ds,  CAPE-Hs  and  CAPE-Ts.  Photos  of  these  vessels  and 
their  principal  characteristics  can  be  found  at  the  results  section  of  Chapter  IIL  Except  for 
the  difference  in  displacement  between  the  different  classes  of  ships,  there  is  another 
significant  difference,  namely  in  the  ramp  design  and  placement  aboard  the  ship  -  stern  or 
side  ramps.  Equipment  is  going  to  be  offloaded  aboard  a  discharge  facility  that  is 
specifically  designed  for  that  purpose.  The  RO/RO  Discharge  Eacility  (RRDE)  is  the 
interface  between  the  ship’s  ramp  and  lighters/landing  craft.  Operational  requirements  for 
equipment  transfer  are  set  at  sea  state  3.  However,  the  current  design  of  the  Ramp/RRDE 
system  is  only  capable  of  operations  in  sea  state  2  and  below.  The  most  limiting 
conditions  result  from  high  relative  motions  between  the  ship  and  the  RRDE  as  well  as 
the  induced  stresses  that  may  exceed  the  yield  stress  of  the  material. 

Successful  implementation  of  a  Ramp/RRDE  Motion  Compensation/Mitigation 
System  will  provide  a  significant  increase  in  military  readiness  and  capability  of  the 
RO/RO  ships  by  enabling  ramp  operations  during  sea  state  3,  thus  increasing  the 
operational  window  by  up  to  50%  (Applebee  et  al.  1997).  Previous  studies  on  the  same 
subject,  done  at  the  Naval  Postgraduate  School  [Ref.  3]  and  the  University  of  Notre 
Dame  [Ref  5],  have  suggested  that  active  or  passive  motion  compensation  can  reduce  the 
motions  and  the  stresses  of  the  ramp.  Moreover,  detailed  finite  element  models  of  the 
various  ship’s  ramps  have  been  created  and  tested  for  stress  and  natural  frequencies 
calculations  [Ref.  6].  In  Ref  [5]  the  hydrodynamic  properties  of  the  ships  were 
approximated  by  a  constant  coefficient  mass-spring-damper  system,  and  active  control 
techniques  were  designed  and  analyzed.  In  Ref.  [3]  several  simplifying  assumptions  were 
made  in  the  development  of  the  dynamic  model  of  the  vessels  and  their  interface.  Ship 
motions  were  approximated  by  the  vertical  plane  equations  only,  and  the  ramp  was 
modeled  as  a  generic  second  order  system. 
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B,  OBJECTIVES 

This  thesis  will  attempt  to  overeome  some  of  the  limitations  of  the  previous 
studies.  It  will  focus  on  creating  a  mathematical  dynamic  model  for  the  ship  and  the 
barge  (RRDF).  These  models  are  going  to  be  created  using  the  latest  calculation  tools  for 
marine  vehicles  so  that  vessel  motions  can  be  calculated  with  the  least  number  of 
simplifying  assumptions.  A  finite  element  model  of  the  ship’s  ramp  is  going  to  be  created 
so  that  it  can  be  used  in  both  the  structural  coupling  of  the  two  vessels  and  in  ramp  stress 
calculations.  Once  the  equations  of  motion  of  the  vessels  in  the  frequency  domain  are 
written,  the  two  vessels  are  going  to  be  coupled  with  the  ramp  by  using  a  specified 
coupling  condition.  This  will  allow  calculation  of  the  new  motion  of  the  ships  along  with 
the  displacement  field  of  the  ramp.  In  this  way  the  response  in  regular  waves  will  be 
calculated.  Finally  the  response  in  a  seaway  of  the  coupled  system  will  be  expressed  in 
terms  of  standard  seaway  parameters. 

For  demonstration  purposes  the  CAPE-D  class  of  ships,  with  stern  ramp,  was 
chosen  along  with  the  new  RRDF  2000  model.  The  coupling  condition  chosen  for 
demonstration  will  be  spring-damper  coupling  condition  with  constant  parameters, 
modeling  a  standard  motion  isolation  system.  The  method  presented  in  this  work  is, 
however,  general  enough  so  that  additional  coupling  conditions  can  be  considered. 
Finally,  a  design  methodology  for  proper  motion  mitigation  is  presented  along  with 
parametric  studies  for  different  isolator  properties. 
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II.  HYDRODYNAMIC  MODEL  OF  THE  SHIP  AND  THE  BARGE 
(RRDF)  RESPONSE  IN  REGULAR  WAVES 

A.  INTRODUCTION 

Modern  seakeeping  eomputations  are  used  in  all  aspeets  of  the  engineering  of  the 
marine  environment.  The  eomputations  are  performed  using  a  variety  of  teehniques-from 
simple  strip  theory  to  extremely  eomplex  fully  nonlinear  unsteady  Reynolds  Averaged 
Navier-Stokes  equations  in  the  time  domain. 

The  major  diffieulties  in  seakeeping  eomputations  are  the  nonlinearities.  These 
are  nonlinearities  assoeiated  with  the  fluid  in  the  form  of  viseosity  and  the  veloeity 
squared  terms  in  the  pressure  equation.  The  free  surfaee  eauses  nonlinear  behavior  due  to 
the  nature  of  the  free-surfaee  boundary  eonditions  and  the  nonlinear  behavior  of  the 
ineident  waves. 

B,  DEVELOPMENT  OF  SEAKEEPING  CALCULATIONS 

For  the  interested  reader  a  detailed  historieal  development  of  seakeeping 
ealeulations  ean  be  found  in  Ogilvie  (1977),  Newman  (1978),  Mauro  (1989)  and  Beek  & 
Reed  (2001)  [Ref.  1].  A  short  historieal  approaeh  to  seakeeping  eomputation  will  be 
presented  here  based  on  Baek  &  Reed  (2001)  [Ref.  1]  where  the  most  reeent  eomputation 
methods  where  evaluated. 

The  effort  of  predieting  the  response  of  a  vessel  in  different  sea  eonditions  has  a 
long  history  starting  with  Fronde’s  original  work  on  rolling  (Froude  1861).  Modern 
eomputations  began  with  two  developments  in  the  1950’s.  The  first  was  the  use  of 
random  proeess  theory  to  determine  the  statisties  of  the  ship  responses  in  a  seaway.  The 
seeond  was  the  development  of  the  linear  ship  motion  theories  to  prediet  the  responses  of 
the  ship  to  regular  waves. 

In  1953  St.  Denis  and  Pierson  in  their  seminal  paper  proposed  a  method  to  prediet 
the  statisties  of  ship  responses  to  a  realistie  seaway.  Using  speetral  methods  developed  in 
other  fields,  they  related  the  speetral  density  of  ship  responses  to  the  input  oeean  wave 
speetrum.  Two  of  their  assumptions  are  eritieal:  1)  the  sea  surfaee  is  an  ergodie,  Gaussian 
random  proeess  with  zero  mean  and  2)  the  ship  can  be  represented  by  a  linear  system. 
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Once  the  probability  density  function  for  a  given  response  is  known,  all  the  desired 
statistics  of  the  response  can  easily  be  determined.  The  linear  system  assumption  allows 
the  spectral  density  of  any  given  response  to  be  found  by  multiplying  the  incident  wave 
spectrum  by  the  square  of  the  Response  Amplitude  Operator  (or  RAO)  of  the  desired 
response.  At  any  single  frequeney,  the  RAO  is  the  amplitude  and  the  phase  of  the  desired 
response  to  regular  ineident  waves  aeting  on  the  vessel  at  the  given  frequency. 

In  the  1950’s  the  development  of  analytic  prediction  techniques  was  started.  The 
thin-ship  approximation  of  Michell  (1898)  was  examined  eritically  by  Peters  and  Stoker 
(1957).  The  first  order  theory  was  rather  trivial  in  that  it  balanced  hydrodynamic  forces 
due  to  the  undisturbed  ineident  wave  pressure  field  (the  Froude-Krylov  exciting  force) 
and  the  hydrostatie  restoring  forees  with  the  ship’s  mass  time  aceeleration  term.  Newman 
(1961)  avoided  the  shorteomings  of  Peters  and  Stoker  by  introducing  refinements  using  a 
systematie  expansion  in  multiple  small  parameters.  Unfortunately,  as  with  thin-ship 
theory,  most  nontrivial  hydrodynamic  effects  are  higher  order  compared  to  the  Froude- 
Krylov  exciting  force  and  the  hydrostatic  restoring  foree. 

At  the  same  time  an  alternative  strip  slender-body  was  also  being  studied.  Korvin- 
Kroukovsky  (1955)  did  the  initial  work.  Using  a  combination  of  slender-body  theory  and 
good  physieal  insight,  they  developed  a  theory  for  heave  and  piteh  that  was  suitable  for 
numerical  computations  for  the  newly  emerging  digital  computers.  Strip  theory  was  the 
first  ship  motion  theory  that  gave  results  with  enough  engineering  accuracy  that  the 
predicted  motions  were  useful  for  design  purposes.  In  the  late  1960’s  more 
comprehensive  strip  theories  were  developed,  most  widely  cited  is  Salvensen,  et  al. 
(1970).  Using  a  combination  of  mathematics  and  judicious  assumptions,  these  researchers 
arrived  at  a  form  of  strip  theory  that  today  is  still  the  most  widely  used  method  for 
seakeeping  eomputations  of  ships. 

Strip  theory  is  a  short  wavelength  theory  and  slender-body  theory  is  a  long 
wavelength  theory.  Attempts  have  been  made  to  bridge  the  gap  and  find  the  theory  that  is 
valid  over  a  wide  frequency  range.  The  interpolation  theory  of  Mauro  (1970)  and  the 
unified  theory  of  Newman  are  typical  examples.  Comparisons  with  experimental  results 
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by  Sclavounos  (1990)  have  indicated  improved  predictions  relative  to  strip  theory 
predictions. 

By  the  late  1970’s  the  Neumann-Kelvin  approach  was  starting  to  be  used.  In  this 
approach  the  body  boundary  condition  is  applied  on  the  mean  position  of  the  exact  body 
surface  and  the  linearized  free-surface  boundary  condition  is  used.  The  traditional 
approach  to  solve  the  Neumann-Kelvin  problem  is  to  use  boundary  integral  methods  in 
which  the  solution  is  formulated  in  terms  of  integrals  of  fundamental  singularities  (source 
and  dipoles)  over  the  surface  surrounding  the  fluid  domain. 

Hess  and  Smith  (1964)  pioneered  boundary  element  methods  for  flow  without  a 
free  surface  (equivalent  to  a  double-body  flow  with  a  rigid  body  surface).  Using  just  a 
source  distribution,  they  subdivide  the  body  surface  into  N  flat  quadrilaterals  over  which 
the  source  strength  was  assumed  constant.  The  flat  quadrilaterals  were  often  called  panels 
and  now  the  term  “panel  methods”  has  come  to  mean  any  solution  technique  in  which  the 
body  surface  (and  possibly  other  surfaces  of  the  problem)  has  been  subdivided.  In  panel 
methods,  two  tasks  require  almost  all  the  computational  effort.  The  first  is  setting  up  the 
influence  matrix  that  requires  multiple  evaluations  of  the  Green  function  for  the  problem. 
The  second  is  solving  the  resulting  linear  system  of  equations. 

Boundary  element  methods,  while  the  most  popular,  are  not  the  only  methods 
available  to  solve  the  Neumann-Kelvin  problem.  Finite  element  or  finite  difference 
approaches  can  also  be  used.  The  only  drawback  is  the  mesh  generation  for  a  complex 
body  geometry.  The  total  computation  effort  and  accuracy  of  the  different  numerical 
techniques  depends  on  the  details  of  the  problem,  the  code,  and  the  computers  on  which  it 
is  to  be  run.  Several  commercial  codes  are  available,  the  first  probably  being  Garrison 
(1978)  and  the  most  widely  used  is  Wave  Analysis  Massachusetts  Institute  of  Technology 
(WAMIT)  (Korsmeyer  et  al.  1988).  The  codes  have  been  extended  to  include  second 
order  mean  drift  and  slowly  vary  forces. 

WAMIT  is  a  radiation/diffraction  program  developed  for  the  analysis  of  the 
interaction  of  surface  waves  with  offshore  structures.  WAMIT  was  used  for  the 
mathematical  modeling  of  the  barge  and  the  ship  and  is  based  on  a  three-dimensional 
panel  method.  The  theory  and  formulation  of  WAMIT  is  described  in  the  following 
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section.  The  most  recent  developments  and  computation  methods  are  described  in  detail 
in  [Ref.  1] 

C.  COMPUTATION  OF  MOTIONS  IN  REGULAR  WAVES 

The  response  of  a  ship  advancing  in  a  seaway  is  a  complicated  phenomenon 
involving  the  interactions  between  the  vessel  dynamics  and  several  distinct 
hydrodynamic  forces.  All  ship  responses  are  nonlinear  to  some  extent  but  in  many  cases 
when  nonlinearities  are  small  a  linear  theory  will  yield  good  predictions.  As  discussed  in 
the  previous  section,  the  seaway  can  be  considered  a  random  process  and  spectral 
techniques  can  be  used  to  define  the  characteristics  of  the  seaway.  By  knowing  the 
responses  of  a  ship  to  regular  waves  of  different  frequencies,  we  can  predict  the  statistics 
of  the  response  to  actual  random  seaways.  In  Ref.  [2]  the  interested  reader  could  find 
more  details  about  the  computation  of  vessel’s  motion,  here  a  short  discussion  will  be 
presented. 

1,  Calculation  Of  Motions  Using  Six  Degrees  Of  Freedom 

A  ship  advancing  at  a  steady  mean  forward  speed  with  arbitrary  heading  in  a  train 
of  regular  waves  will  move  in  six  degrees  of  freedom.  That  is  the  ship’s  motion  can  be 
considered  to  be  made  up  of  three  translational  components,  surge,  sway  and  heave,  and 
three  rotational  components,  roll,  pitch  and  yaw  as  shown  in  Figure  1.  The  six  modes  of 
motion  are  therefore: 

=  surge  =  translation  in  x-direction 

^2  =  sway  =  translation  in  y-direction 

^3  =  heave  =  translation  in  z-direction 

^4  =  roll  =  rotation  about  x-axis 

^5  =  pitch  =  rotation  about  y-axis 

=  yaw  =  rotation  about  z-axis 
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In  the  problem  to  be  formulated  here,  as  origin  will  be  considered  a  point  in  the 
waterline  (z  =  0),  on  the  centerline  (j  =  0)  and  at  the  longitudinal  center  of  gravity 
(which  is  also  the  longitudinal  center  of  buoyancy  for  a  freely  floating  body). 

In  order  to  be  able  to  compute  all  responses  of  any  vessel  to  regular  waves,  it  is 
necessary  to  deal  with  the  full  six  degrees  of  freedom  equations,  considering  important 
couplings  among  them.  Here  only  the  final  linear  form  of  the  equations  of  motion  will  be 
presented.  Detailed  derivations  of  the  equations  of  motion  may  be  found  in  Salvensen,  et 
al  (1970),  Newman  (1977),  Ogilvie  (1964),  or  Wehausen  (1971). 

The  linearization  of  the  equations  is  made  on  the  basis  of  small  motions.  The 
motions  will  in  general  be  small  if  the  ship  is  stable  and  the  incident  wave  amplitude  is 
relatively  small.  The  principal  exceptions  to  this  rule  are  resonant  situations  where 
damping  is  small,  e.g.,  roll  resonance  in  beam  seas.  Experimental  and  theoretical 
investigations  have  shown  that  a  linear  analysis  of  ship  motions  gives  excellent 
predictions  over  a  wide  variety  of  sea  conditions  and  vessel  types.  In  situations  where  the 
linear  theory  assumptions  are  in  doubt,  only  experiments  can  be  used  to  compare  the 
calculation  results. 

In  linear  theory,  the  responses  of  the  vessel  will  be  linear  with  (i.e.  directly 
proportional  to)  wave  amplitude  and  occur  at  the  frequency  at  which  the  ship  perceives 
the  incident  waves.  This  is  an  important  point  because  even  though  it  allows  the  use  of 


powerful  linear  tools,  it  also  puts  limitations  on  the  results,  whieh  must  be  reeognized. 
Sinee  only  the  response  of  the  vessel  in  sinusoidal  waves  (regular  waves)  is  being 
eonsidered  in  this  seetion,  the  vessel’s  responses  in  the  time  domain,  will  be 

sinusoidal  at  the  frequeney  of  eneounter  and  ean  be  written  as  follows: 

=  7=1,2...6  (1) 

where,  is  the  eomplex  amplitude  of  the  vessel’s  response  in  the  y-th  direetion  and 
7  =1,2,3,. ..,6  refer  to  surge,  sway,  heave,  roll,  piteh,  and  yaw  respeetively  and  co^  is  the 
frequeney  of  eneounter  and  equals 

2 

cOq - -U^co^jii  (2) 

g 

where,  y/  is  the  angle  that  the  ship  moves  with  respeet  to  the  regular  waves,  co^  is  the 
frequeney  of  the  regular  wave  and  f/gis  the  speed  of  the  ship.  When  the  ship’s  speed  is 
zero,  whieh  is  the  ease  eonsidered  here,  then  co^  equals  co^ . 

Using  the  linear  theory,  the  motion  of  a  vessel  ean  be  modeled  as  a  generalized 
spring-mass-damper  system.  The  equations  of  motion  then  are  written  as  six 
simultaneous  linear  equations  in  the  time  domain: 

=  y  =  l,2,...6  (3) 

k=\  k=\  k=\ 


Using  the  linear  theory  assumption  and  Equation  (1),  f}^  and  %  can  be  expressed  as: 


(4) 


Substituting  Equations  (4)  into  Equation  (3)  and  eliminating  the  e'®"'  term  we  have  the 
system  of  equations  in  the  frequency  domain: 


I[- 

k=\ 


co: 


)  + 


(5) 
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The  Aj^.  terms  are  ealled  hydrodynamic  added  mass.  These  terms  are  in  phase  with 
vertieal  aeeelerations  and  physieally  Aji^  represent  the  foree  eomponent  in  the  j  -th  mode 
of  motion  due  to  the  aeeeleration  in  the  k  -th  mode  of  motion.  The  terms  eorrespond 
to  hydrodynamic  damping  in  phase  with  vertieal  veloeity,  and  physieally  Bj^.  represent 
the  foree  eomponent  in  the  j  -th  mode  of  motion  due  to  the  veloeity  in  the  k  -th  mode  of 
motion.  Generally  both  A^^.  and  are  funetions  of  frequeney.  Terms  involving  the 
eoeffieients  are  the  restoring  forces  and  moments  representing  the  net  hydrostatie 
buoyaney  effeets  of  the  ship  motions  and  similarly  Cj^.  are  the  hydrostatie  restoring  foree 
eoeffieients  in  the  j  -th  direetion  due  to  k  -th  motion.  The  terms  represent  the 

hydrodynamic  wave  exciting  forces  and  moments  and  are  usually  subdivided  into  two 
eomponents  as  follows: 

F^=F‘+F^  (6) 

The  terms  Fj  are  the  eomplex  amplitudes  of  the  wave  exeiting  forees  and  moments  due 

to  incident  waves,  also  known  as  Froude-Krylov  exciting  forces  after  the  classical  work 
on  ship  rolling  done  by  Froude  (1861),  and  generalized  to  six  degrees  of  freedom  by 
Krylov  (1898).  The  Froude-Krylov  excitations  represent  the  integration  of  the  pressure 
over  the  body  surface  that  would  exist  in  the  incident  wave  system  if  the  body  were  not 
present.  The  F^  terms  are  the  complex  amplitudes  of  the  wave  exciting  forces  and 

moments  due  to  diffracted  waves.  The  diffraction  exciting  forces  and  moments  are 
caused  by  the  diffraction  or  modification  of  the  incident  waves  due  to  the  presence  of  the 
vessel.  The  computation  of  these  terms  is  possible  by  using  the  so-called  Haskind 
Relations,  which  are  presented  later  in  this  Chapter.  The  Froude-Krylov  forces  and 
moments  are  sometimes  used  to  approximate  the  total  exciting  forces  since  this  is  a 
considerably  simpler  task.  This  is  a  good  approximation  only  if  the  wavelength  is  much 
larger  than  the  length  of  the  vessel,  such  as  in  the  case  of  a  small  underwater  vehicle  in 
waves.  For  shorter  wavelengths  the  approximation  is  increasingly  inaccurate  because  the 
diffraction  forces  become  significant. 
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Figure  2.  Sketch  Defining  the  Relative  Position  and  Geometry  of  the  CAPE-D  and 

RRDF 


Equations  (5)  hold  for  a  vessel  floating  in  the  open  sea  without  any  object  in  the 
proximity.  In  the  case  under  consideration  here,  a  ship  (CAPE-D  or  CAPE-H)  and  a 
barge  (RRDE)  are  floating  very  close  to  one  another  while  maintaining  zero  speed.  The 
relative  position  and  the  geometry  of  the  two  vessels  can  be  seen  in  Figure  (2).  The  short 
distance  between  the  two  vessels  stimulates  radiative  effects.  This  is  the  result  of  the 
waves  striking  on  one  vessel,  bouncing  off  its  hull,  and  then  striking  the  other  vessel. 
Considerable  coupling  is  created  between  the  two  vessels,  which  depends  on  the  relative 
direction  of  the  incident  waves.  In  order  to  model  accurately  the  motion  of  the  ship  and 
the  barge,  coupling  coefficients  have  to  be  introduced,  and  the  equations  of  motion  of  the 
two  vessels  have  to  be  solved  simultaneously.  Now  the  system  is  made  of  twelve 
equations  and  twelve  unknowns.  The  twelve  unknowns  are  the  six  degrees  of  freedom  of 
the  ship  and  the  six  degrees  of  freedom  of  the  barge.  In  matrix  form  the  system  can  be 
represented  by  a  single  matrix  equation 

(-a>l[M  +  A]  +  il0„  [ii]  +  [C]){^}  =  {F)  (7) 

From  Equation  (7)  and  by  knowing  the  values  of  the  coefficients  and  the  exciting 
forces,  we  can  obtain  the  twelve  modes  of  motion  for  every  frequency  as  follows: 
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{n]  =  (-OJ>l[M  +  A\  +  iQ),[B]  +  [C])  '{F} 


(8) 


The  solution  of  the  above  equation  provides  the  RAO  for  every  motion  of  the 
ship’s  and  the  barge’s  center  of  gravity.  Knowing  the  motions  of  the  gravitational  center 
and  by  observing  the  geometry,  the  vertical  motion  of  every  point  aboard  the  ship  or  the 
barge  can  be  found  as: 

^  =  %  +  (9) 

where  ^  is  the  complex  local  amplitude  of  the  vertical  motion  of  the  given  point  and  x,  y 
are  its  coordinates  on  the  body  fixed  frame. 

The  determination  of  the  coefficients,  the  amplitudes  of  the  exciting  forces  and 
moments  in  Equation  (8)  represents  the  major  task  in  ship  motions  calculations.  The 
problem  can  be  simplified  by  the  use  of  advanced  computational  techniques  as  those 
discussed  in  the  previous  section. 

2.  Calculation  Of  Motions  Using  a  Three  Degrees  Of  Freedom 
Approximation 

Many  times  in  calculating  the  motions  of  a  vessel  can  be  assumed  that  the 
coupling  between  the  six  degrees  of  freedom  is  not  strong  and  Equations  (5)  can  be 
decoupled.  This  is  valid  for  the  case  of  an  unrestrained  ship  with  port/starboard  symmetry 
where  the  six  equations  may  be  uncoupled  into  two  sets  of  three  equations.  The  vertical- 
plane  or  longitudinal  motions  (surge,  heave  and  pitch)  are  uncoupled  from  the  horizontal- 
plane  or  transverse  motions  (sway,  roll,  and  yaw).  It  should  be  noted  that  the  lack  of 
coupling  between  the  vertical  and  horizontal  modes  is  a  consequence  of  linear  theory.  In 
nonlinear  theories  such  cross-coupling  may  be  present.  Eor  some  ship  motion  problems 
this  non-linear  coupling  can  be  very  important.  Eor  example,  there  is  a  non-linear  heave- 
roll  cross-coupling  that  can  lead  to  roll  instabilities  and  eventual  ship  capsizing  (Kerwin, 
1955)  or  (Ogilvie  and  Beck,  1973).  Another  example  is  the  nonlinear  pitch-yaw  coupling 
that  results  from  varying  submergence  of  the  bow  due  to  pitch  motion  (Korvin- 
Kroukovsky,  1980). 
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In  Ref.  [3]  a  detailed  formulation  of  the  equations  (using  the  above 
approximation)  of  the  ship  and  the  barge,  under  eonsideration,  can  be  found.  Proximity 
effects  were  included  into  the  equations,  thereby  introducing  hydrodynamic  coupling 
between  the  barge  and  the  ship.  The  detailed  formulation  will  not  presented  here, 
however  some  of  the  results  obtained  by  this  method  are  going  to  be  presented  in  the 
results  section  for  comparison. 

3,  Determination  Of  the  Hydrodynamic  Forces  and  Coefficients  Using 
WAMIT 

In  the  present  study  we  utilized  results  obtained  from  WAMIT.  A  short  discussion 
of  how  WAMIT  calculates  the  coefficients  under  determination  will  be  presented  here. 
For  more  details  the  reader  should  refer  to  WAMIT  user  manual  [Ref.  4]. 

The  objective  of  WAMIT  is  to  evaluate  the  unsteady  hydrodynamic  pressure, 
loads  and  motions  of  the  body,  as  well  as  the  induced  pressure  and  velocity  in  the  fluid 
domain.  The  free-surface  and  body-boundary  conditions  are  linearized,  the  flow  is 
assumed  to  be  potential,  free  of  separation  or  lifting  effects.  A  harmonic  time  dependence 
is  adopted.  Figure  3  illustrates  a  three-dimensional  body  interacting  with  plane 
progressive  waves  in  water  of  finite  water  depth  H 


Figure  3.  Sketch  Defining  the  Coordinates,  Fluid  Depth,  and  Wave-Heading  Angle 

from  Ref.  [3]. 
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The  Cartesian  eoordinate  system  (x,y,z)  defined  in  Figure  3  is  stationary  relative 
to  the  undisturbed  position  of  the  free  surfaee  and  body.  Its  origin  may  lie  on,  above  or 
below  the  free  surfaee,  as  the  partieular  applieation  requires.  The  body  geometry  input  to 
WAMIT  is  defined  relative  to  that  system.  Here,  the  origin  is  loeated  on  the  free  surfaee 
for  the  sake  of  simplieity. 

The  assumption  of  a  potential  flow  permits  the  definition  of  the  flow  veloeity  as 
the  gradient  of  the  veloeity  potential  satisfying  the  Laplaee  equation 

¥"^>  =  0  (10) 

in  the  fluid  domain.  The  harmonie  time  dependenee  allows  the  definition  of  a  eomplex 
veloeity  potential  ^ ,  related  to  by 

^  =  Re(^e“^‘)  (11) 


where  Re  denotes  the  real  part,  co^  is  the  frequeney  of  the  ineident  wave  and  t  is  time. 

The  ensuing  boundary-value  problem  will  be  expressed  in  terms  of  the  eomplex  veloeity 
potential,  <p  with  the  understanding  that  the  produet  of  all  eomplex  quantities  with  the 

faetor  applies.  The  linearized  form  of  the  free-surfaee  eondition  is 

(p^-K(p  =  0,  z  =  0  (12) 


where  K  =  of  !  g  ,  and  g  is  the  aeeeleration  of  gravity.  The  veloeity  potential  of  the 
ineident  wave  is  defined  by 


_  igA  cosh[v(z  +  //)] 

°  CO  cosh  vH 

where  the  wavenumber  v  is  the  real  root  of  the  dispersion  relation 


(13) 


—  =  vtanhv//  (14) 


and  (5  is  the  angle  between  the  direction  of  propagation  of  the  incident  wave  and  the 
positive  X  axis  as  defined  in  Figure  3. 
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The  linearization  of  the  problem  permits  the  deeomposition  of  the  veloeity 
potential  into  the  radiation  and  diffraetion  eomponents 


(P=(PR+(PD 

(15) 

6 

(p„=io)Y,Vj<Pj 

(16) 

(Pd  =  (Po  +  (Pi 

(17) 

The  eonstants  rjj  denote  the  eomplex  amplitudes  of  the  body  oseillatory  motion 
in  its  six  rigid-body  degrees  of  freedom,  and  (Pj  the  eorresponding  unit-amplitude 
radiation  potentials.  The  veloeity  potential  (p^  represents  the  scattered  disturbanee  of  the 
ineident  wave  by  the  body  fixed  at  its  undisturbed  position.  The  sum  (17)  will  refer  to  as 
the  diffraction  potential  tpj^ . 

On  the  undisturbed  position  of  the  body  boundary,  the  radiation  and  diffraetion 
potentials  are  subjeet  to  the  eonditions 


(p.  =n 

(18) 

o 

II 

s: 

(19) 

where  [n^,n2,nf\  =  n  and  [n^,n^,nf^=x'><n,  x  =  [x,y,z).  The  unit  veetor  n  is  normal  to 
the  body  boundary  and  points  out  of  the  fluid  domain. 

In  WAMIT  the  boundary  value  problems  (10)-(19)  are  solved  by  using  Green’s 
theorem  to  derive  integral  equations  for  the  radiation  and  diffraetion  veloeity  potentials 
on  the  body  boundary.  The  integral  equation  satisfied  by  the  radiation  veloeity  potentials 
tp  -  on  the  body  boundary  takes  the  form 

iTitPi  (x)  +  (p.  np  {<^;x)d<^  (20) 

where  denotes  body  wetted  surfaee  at  ealm  water  and  G{xp)  is  the  Green  function. 


The  corresponding  equation  for  the  total  diffraction  velocity  potential  tp^  is 
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iTTCPo  +  =  4m  (^)  (21) 

The  diffraction  potential  may  also  be  obtained  from  Equation  (17)  after  solving 
for  the  scattered  potential  cp^ .  The  equation  for  the  scattered  velocity  potential  is 

=  (22) 

The  Green  function  G{x;^)  is  referred  to  as  the  wave  source  potential.  It  is  the 
velocity  potential  at  the  point  x  due  to  a  point  source  of  strength  -4;r  located  at  point  ^ . 
It  satisfies  the  free-surface  and  radiation  conditions,  and  in  infinite  water  depth  is  defined 
by 

G(x-Q  =  -  +  -  +  —\dk——J,(kR)  (23) 

r  r  ;r  k-K 

=  +  {y-rif  +  {z-Cf  (24) 

r'^={x-^f  +  [y-ilf  +  {z  +  Cf  (25) 


where  (x)  is  the  Bessel  function  of  zero  order.  In  finite  depth,  the  Green  function  is 
defined  by 


I  I  c~  (A:  +  7f)coshA:(z  +  7/)coshA:(^  +  7/) 

— +  — +  2  dk - e 

r  r  •'0  A:  sinh  cosh 


J,{kR) 


=  {x-^f  +  {z  +  C  +  2Hf  (27) 


In  both  expressions  (23)  and  (26)  the  Fourier  k  -integration  is  intended  above  the 
pole  on  the  real  axis  in  order  to  enforce  the  radiation  condition.  Efficient  algorithms  for 
the  evaluation  of  the  infinite  and  finite-depth  wave-source  potentials  and  their  spatial 
derivatives  have  been  developed.  Provision  has  been  made  in  WAMIT  to  permit  the 
logarithmic  singularity  and  its  derivatives  to  be  integrated  analytically  in  the  solution  of 
the  integral  equations  when  the  source  and  field  points  are  close  to  each  other  and  to  the 
free  surface. 
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As  already  mentioned  in  the  above  section,  WAMIT  utilizes  the  panel  method  for 
discretization.  This  allows  for  numerical  solution  of  the  above  integral  equations.  The 
mean  position  of  the  body’s  wetted  surface  is  approximated  by  a  collection  of 
quadrilaterals.  Four  vertices,  lying  on  the  body’s  surface,  define  each  quadrilateral.  In 
general  the  quadrilaterals  defined  above  are  not  plane,  but  if  a  sufficiently  fine 
discretization  is  used  for  the  boundary  surface  with  continuous  curvature,  each  element 
will  approach  a  plane  surface. 

The  above  solution  of  the  velocity  and  scattering  potentials  allows  WAMIT  to 
calculate  the  coefficients  of  Equation  (5),  which  can  be  written  in  matrix  form  as 

(-0)^,  [M  +  A]  +  ico,  [B]  +  [C]){fj}  =  {F}  (28) 


where  M  is  the  inertia  matrix  defined  as 


m 

0 

0 

0 

mz^ 

-my^ 

0 

m 

0 

-mz^ 

0 

mx^ 

0 

0 

m 

my^ 

-mx^ 

0 

0 

-mz^ 

my^ 

In 

In 

In 

mz^ 

0 

-mx^ 

Ill 

hi 

-my^ 

mx^ 

0 

hi 

hi 

(29) 


m  is  the  body  mass,  {xg,yg,^g)  are  the  coordinates  of  the  center  of  gravity  and  I^j  are 
the  moments  of  inertia  in  terms  of  the  corresponding  radii  of  gyration. 


The  forces  and  other  quantities  evaluated  by  WAMIT  are  output  in  a  standard 
nondimensional  form,  in  terms  of  the  appropriate  combinations  of  the  water  density  p , 
the  acceleration  of  gravity  g ,  the  incident-wave  amplitude  A  ,  the  frequency  CO ,  and  the 
length  scale  L . 


All  hydrostatic  data  can  be  expressed  in  the  form  of  surface  integrals  over  the 
mean  body  wetted  surface  5*^ ,  by  virtue  of  Gauss’  divergence  theorem.  So  the  mass  can 
be  found  by 


m  =  p\f 


where  V  is  the  underwater  volume  and  can  be  calculated  by 

16 


(30) 


(31) 


V  =  n^xdS  n2ydS  =  ~||^  n^zdS 

The  elements  of  the  hydrostatic  and  gravitational  eoefficients  matrix  C  in 
Equation  (28)  ean  be  found  as  follows: 

C(3,3)  =  /7gJ£^V5 

C[X^)  =  pgW^yn^dS 

C(3,5)  =  -/7gjj^^  xn^dS 

C(4, 4)  =  /7g y^n^dS  +  pg^z^  - mgz^ 

C(4,5)  =  -/7gj£^xjn3J5 
C  ( 4, 6 )  =  -pg^Xi,  +  mgXg 
C(5,5)  =  /7gJf  x^n2dS  +  pg^Zf^-rngz^ 

C[5,6)  =  -pg\/y,  +  mgy^ 


WAMIT’s  output  format  is  in  non-dimensional  form  that  is  defined  as: 

C(3,3)  =  C(3,3)//7gE^ 

C(3,4)  =  C(3,4)//7gE^ 

C(3,5)  =  C(3,5)//7gZ^ 

C(4,4)  =  C(4,4)//7gE^ 

C[A,5)  =  C[A,5)I  pgL^ 

C(4,6)  =  C(4,6)//7gE" 

C{5,5)  =  C[5,5)I  pgL^ 

C{5,6)  =  C{5,6)I  pgL^ 


where  C{i,j)  =  C[j,i)  for  all  i,j  exeept  for  C(4,6)  and  C (5, 6) .  For  all  other  values 
of  the  indices  /,7,C(/,7)  =  0.  In  particular,  C(6,4)  =  C(6,5)  =  0. 

The  added  mass  and  added  damping  coefficients  of  the  matrices  A  and  B  of 
Equation  (28)  can  be  found  by  the  following  equations: 
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(34) 


in  the  non-dimensional  output  form 


(35) 


where  k  =  3  for  /,7  =  1,2,3,  k  =  A  for  /  =  1,2, 3,7  =  4,5,6  or  /  =  4, 5,6,7  =  1,2,3 
k  =  5  for  /,7  =  4,5,6. 


The  exeiting  forees  are  ealeulated  both  from  the  Haskind  Relations 


X,=-iO)p\\[n.(p,-(p, 


dn  j 


dS 


(36) 


and  from  direet  integration  of  the  hydrodynamie  pressure 

X.  =  -iO}p\\^  n,(podS 


(37) 


The  non-dimensional  output  format  is  defined  as: 


X,.  = 


X. 


pgAVr 


(38) 


where  m  =  l  for  /  =  1, 2, 3  and  m  =  3  for  z  =  4, 5, 6 . 

D,  RESULTS 

Numerieal  results  are  presented  for  the  Cape-D  elass  of  ships  and  the  RRDF2000 
(Roll-On  Roll-Off  Diseharge  Faeility).  The  RRDF  is  made  up  of  a  number  of 
intereonneeted  modules,  in  an  asymmetrieal  way,  as  ean  be  seen  in  Figure  4.  Eaeh 
module  has  8x8x40  ft  nominal  dimensions,  and  floats  at  an  approximate  draft  of  2  ft.  The 
RRDF  is  plaeed  astern  of  the  Cape-D  ship  whieh  is  equipped  with  a  stem  ramp  only.  The 
Cape-D  ship  has  a  length  of  634  ft,  beam  97  ft,  nominal  draft  32.4  ft,  and  displaees 
approximately  36,000  tons.  This  ship  is  part  of  the  Ready  Reserve  Foree  and  is  one  of  the 
main  ships  used  for  dry  eargo  transfer  by  the  Military  Sealift  Command’s  (MSC) 
Strategie  Sealift  Forees.  A  pieture  of  the  ship  is  shown  in  Figure  5.  Raw  WAMIT  results 
for  this  eonfiguration  in  terms  of  added  mass,  damping,  hydrostatie,  and  wave  exeiting 

forees  were  provided  by  NSWC,  Carderoek  Division.  Caleulations  were  performed  for  30 
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frequencies  ranging  from  0.3  to  2.5  rad/sec  and  for  24  headings,  0  to  345  degrees  in  15 
degree  increments.  We  remind  the  reader  that,  as  shown  in  Figure  3,  0  deg  corresponds  to 
following  (astern)  seas,  180  deg  to  head  seas,  90  deg  to  starboard  beam  seas,  and  270  deg 
to  port  beam  seas. 


Figure  5.  Typical  Views  of  the  Cape-D  Ship  Class. 


19 


Regular  wave  results  are  presented  in  Figures  A1  through  A44  in  Appendix  A. 
Figures  A1  and  A2  show  the  heave  response  for  the  Cape-D  per  unit  wave  amplitude  for 
various  wave  direetions.  It  should  be  emphasized  that,  at  this  point,  the  ramp  is  not 
ineluded  in  the  formulation.  Therefore,  there  is  no  struetural  coupling  between  the  two 
bodies,  although  of  course  hydrodynamic  coupling  is  fully  included.  The  results  are 
port/starboard  symmetric.  This  shows  that  ship  response  is  unaffected  by  the  presence  of 
the  RRDF,  a  result  which  is  consistent  with  physical  intuition.  Figures  A3  and  A4  show 
the  corresponding  heave  response  for  the  RRDF.  It  can  be  seen  that  in  this  case,  a  slight 
asymmetry  between  port  beam  and  starboard  beam  seas  has  developed.  This  is  due  to  the 
asymmetric  nature  of  the  RRDF  as  shown  in  Figure  4.  Also,  heave  response  is  less  for 
head  seas  than  following  seas,  a  result  which  is  partially  due  to  the  sheltering  of  the 
RRDF  by  the  ship.  As  expected  (Ref  [8])  the  roll  motion  for  a  large  monohull,  as  the 
CAPE-D,  is  lightly  damped  especially  for  low  frequencies.  This  explains  the  large  values 
for  the  ship’s  roll  motion  at  wave  periods  greater  than  18  seconds.  Corresponding  plots, 
with  similar  conclusions,  for  the  roll  and  pitch  responses  are  shown  in  Figures  A5 
through  A12.  Only  heave,  pitch,  and  roll  motions  are  presented  since  these  motions 
determine  the  vertical  motion  at  the  interface  between  ship  and  RRDF,  which  is  analyzed 
in  the  following  chapters. 

Figures  A13  through  A28  are  presented  in  order  to  gain  some  understanding  on 
the  hydrodynamic  coupling  aspects  of  the  motions  of  the  two  bodies.  These  figures  show 
a  comparison  between  the  full  12  degrees  of  freedom  (6  from  each  body)  solutions  and 
simplified  reduced  approximations  where  only  vertical  plane  motions  are  considered.  The 
latter  approximation  is  often  utilized  in  seakeeping  studies,  since  it  is  generally  accepted 
that  linear  vertical  plane  motions  (heave,  pitch,  and  roll)  are  decoupled  from  horizontal 
plane  motions  (sway,  yaw,  and  surge).  Both  the  amplitude  and  the  phase  angle  are 
included  in  the  figures.  Based  on  this  set  of  figures,  the  following  conclusions  can  be 
made: 

1.  Ship  vertical/horizontal  plane  coupling  is  negligible  for  following  seas,  and  it  is 
small  for  head  seas. 

2.  Ship  vertical/horizontal  plane  coupling  can  be  significant  for  beam  seas. 
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3.  RRDF  coupling  between  horizontal  and  vertical  plane  motions  is  significant 
throughout  the  wave  direetion  and  frequeney  range  and  eannot  be  negleeted.  This 
means  that  any  attempts  to  simplify  the  motion  of  the  RRDF  into  an  equivalent 
low  order  dynamieal  system  will  have  to  proeeed  very  eautiously.  In  this  thesis, 
we  will  continue  to  utilize  the  full  12-degree  of  freedom  hydrodynamic  model. 

Figures  A29  through  A44  present  the  vertical  motion  at  the  two  eonnection  points 
between  the  ship  and  the  RRDF.  Point  A  is  the  starboard  eonneetion  point,  and  point  B  is 
at  port.  These  motions  are  used  as  the  basis  for  the  struetural  eoupling  in  the 
Ship/ramp/RRDF  system  and  also  form  some  of  the  primary  Response  Amplitude 
Operators  (RAO)  that  are  further  studied  in  the  next  chapters.  The  phase  angle  for  the 
Ship/RRDF  vertical  motions  is  also  shown  in  the  figures,  in  order  to  explain  the  fact  that 
the  relative  motion  between  the  two  points  is  not  simply  the  difference  between  the 
amplitudes  of  the  two  vertical  motions.  Points  where  the  relative  motion  experiences  a 
peak  will  become  very  important  in  the  subsequent  chapters,  since  they  will  relate  to 
excessive  twist  and  stresses  in  the  ramp. 
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III.  RAMP  MODELING 


A,  INTRODUCTION 

Following  the  numerical  solution  for  the  motions  of  the  ship  and  the  barge,  the 
formulation  of  the  mathematical  model  of  the  ramp  between  the  ship  and  the  barge  must 
be  done.  Previous  work  done  on  this  subject  has  shown  shortcomings  that  have  to  be 
overcome.  In  Ref.  [3]  a  simplified  second  order  system  was  used  to  simulate  the  ramp, 
also  the  solution  of  the  barge  and  the  ship  motion  was  found  using  the  three  degrees  of 
freedom  hydrodynamic  approximation  described  in  the  previous  chapter.  In  Ref.  [5],  an 
analytical  model  for  the  ramp  was  created  by  using  Timoshenko’s  plate  theory  but  a 
simplified  spring-mass-damper  system  was  used  for  the  barge.  In  Ref  [6]  a  detailed  finite 
element  model  of  the  ramp  was  used  to  show  the  influence  of  the  roll  motion  in  the  stress 
profile  of  the  ramp,  but  the  plate  was  not  directly  coupled  with  the  barge  and  the  ship. 

B,  RAMP  MODELING  USING  FINITE  ELEMENT  METHOD 

In  this  work,  we  will  approximate  the  ramp  by  a  rectangular  plate  with  uniform 
properties.  This  approximation  is  valid  here  since  our  primary  concern  is  to  establish  the 
level  of  influence  that  a  structural  component  such  as  a  ramp  has  on  a  floating  body  such 
as  the  RRDF.  This  approximation  allows  us  to  use  relatively  simple  algorithms  and 
computer  programs  that  can  be  easily  adjusted  in  order  to  examine  different  coupling 
conditions.  For  comparison  purposes  only  the  relative  magnitudes  of  the  stresses  are  of 
interest  for  the  different  coupling  conditions,  and  not  their  exact  values.  Once  the 
ramp/RRDF  influence  levels  have  been  established,  one  could  use  more  detailed  finite 
element  techniques  such  as  those  presented  in  [6]  for  accurate  stress  calculations,  or  the 
strategies  developed  in  [5]  for  active  control  studies. 

A  finite  element  model  for  the  ramp  was  created  using  shear  deformable  plate 
elements,  as  described  in  Ref.  [7].  The  Mindlin/Reissner  theory  was  utilized  for  the 
formulation  of  the  stiffness  matrix.  This  theory  includes  the  effect  of  the  transverse  shear 
deformation  like  the  Timoshenko  beam  theory.  Hence,  a  plane  normal  to  the  midplane  of 
the  plate  before  the  deformation  does  not  remain  normal  to  the  midplane  any  longer  after 
the  deformation.  An  expression  for  the  internal  energy  of  the  shear  deformable  plate 
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should  include  transverse  shear  energy  and  bending  energy  as  well.  The  internal  energy  is 
expressed  as: 


where 

K}=K 

k}=K  r.yY 

are  the  bending  stress  and  strain  eomponents  while, 

K}=K 

k}={nz  Ty^y 

are  the  transverse  shear  components.  Additionally,  /r  is  the  shear  energy 
factor  equal  to  ^  and  V  indieates  integration  over  the  volume. 

Substitution  of  the  equations  (40)-(43)  into  (39)  yields 

U  =  \l{e,Y[D,]{e,]dV+'^l{e,Y[D,\{e,]dV 

in  which 
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E 

1-v" 
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is  the  constitutive  equation  for  the  plane  stress  eondition  and 


[A] 


G  0 
0  G 


Further,  G  is  the  shear  modulus  and  V  is  the  Poisson’s  Ratio. 


(39) 

(40) 

(41) 

(42) 

(43) 
correction 

(44) 

(45) 

(46) 
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Next,  in  order  to  derive  the  element  stiffness  matrix  for  the  shear  deformable 
plate,  the  strains  need  to  be  expressed  in  terms  of  nodal  variables.  The  displacements 
along  the  x,  y  plane  are  given  as 


(47) 

(at) 

(48) 

and  the  transverse  displacement  along  the  z  axis  is 

w  =  w{x,y)  (49) 


where  0^  and  0^  are  the  rotations  of  the  midplane  about  the  y  and  x  axes  respectively. 
The  midplane  is  assumed  to  have  no  inplane  deformation.  For  the  shear  deformable  plate. 


_  dw 
0  = - r 

X  /  XI 

OX 


0  =^-r 

ay 


(50) 

(51) 


where  y  is  the  angle  caused  by  the  transverse  shear  deformation. 


The  displacement  w  and  the  slope  0  are  independent.  Because  of  their 
independence  we  need  shape  functions  to  interpolate  them  independently.  As  a  result,  the 
shear  deformable  plate  element  requires  C°  compatibility.  The  isoparametric  shape 
functions,  ,  are  used  for  the  plate  element  formulation  (Ref.  [7]) 


where 


(52) 

(53) 

(54) 
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(55) 


and  ^,7]  are  the  natural  eoordinates  of  the  isoparametrie  element.  A  point  (^,7)  within 
the  natural  element  is  mapped  into  a  point  {x,y)  within  the  physieal  element,  using  the 
shape  function  described  above,  as 

x=f,H,((,n)x,  (56) 

/■-I 

.’’  =  Z  ".(#.’?)  5'.  (57) 

Then  the  transverse  displacement  and  the  slopes  are  interpolated  as 

w  =  (58) 

=  (59) 

/■-I 

(69) 

i-1 

Here  n  is  the  number  of  nodes  per  element  and  the  same  shape  functions  are  used  for  the 
displacement  and  slope  interpolations.  In  the  formulation  presented  here,  bilinear 
isoparametric  shape  functions  are  used  for  simplicity.  For  more  information  about 
bilinear  isoparametric  elements  and  shape  functions  the  interested  reader  should  read  the 
corresponding  chapter  in  Ref.  [7].  The  bilinear  isoparametric  elements  have  four  nodes 
per  element.  So  the  total  degrees  of  freedom  per  element  are  twelve. 

Both  bending  and  shear  strains  for  each  element  are  computed  from 
displacements 

K)=-z[S,]{rf‘}  (61) 

{^.)  =  lS,]h’}  («) 


where 
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and 


{d-\ 


{(»J,  M,  (»J,  («,),  «'2  "-j  {cy  M,  "'4}' (65) 


and 

Finally,  substitution  of  Equations  (61)  and  (62)  into  the  energy  expression  (44) 
yields  for  eaeh  plate  element 

U  = 

jjs.f  [B,][)j,]dzrfn{rf'}(66) 

As  a  result,  the  element  stiffness  matrix  for  the  plate  bending  ean  be  expressed  as 

where  t  is  the  plate  thiekness  and  J  is  the  Jaeobian  matrix  to  transform  from  the  natural 
eoordinates  (^,7)  to  physieal  eoordinates  (x,y). 

The  mass  matrix  for  eaeh  plate  element  ean  be  found  (Ref.  [7])  as 
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(68) 


[M-]  =  \^,pA\H]'[H]\j\da 

where  /  is  the  length  of  the  individual  element,  p  is  the  mass  density  of  the  ramp,  A  is 
the  eross  seetional  area  of  the  ramp  and  [if]  is  the  diagonal  matrix  of  the  shape 

funetions.  The  above  matrix  is  ealled  the  eonsistent  mass  matrix.  A  more  eomputational 
effieient  form  of  the  above  mass  matrix  is  the  lumped  mass  matrix,  which  is  a  diagonal 
matrix  and  can  be  developed  using  a  method  described  in  Ref  [7].  In  the  case  of  the 
shear  deformable  element  the  lumped  mass  matrix  is  as  follows 

"0  00000000000 

000000000000 

001000000000 

000000000000 

000000000000 

1  000001000000 
M  \  =  ^— 

J  4000000000000 
000000000000 
000000001000 
000000000000 
000000000000 
0  00000000001 

Using  a  program  in  MATLAB  the  mass  matrix  and  the  stiffness  matrix  of  the 
whole  ramp  can  be  found  for  every  boundary  condition  that  the  user  specifies.  The  mesh 
of  the  finite  element  model  consisted  of  100  elements  as  can  be  seen  in  Figure  (6).  A  total 
of  121  nodes  were  used  with  a  total  of  363  degrees  of  freedom. 

C.  VALIDATION  OF  THE  RAMP  MODEL 

The  finite  element  model  of  the  ramp  described  above  approximates  the 
complicated  structure  of  the  physical  ramp,  with  a  plate  of  uniform  properties.  The 
thickness  and  the  modulus  of  elasticity  of  the  model  should  be  selected  so  that  the  natural 
frequencies  of  the  model  are  comparable  with  the  natural  frequencies  of  the  physical 
ramp.  In  Ref  [6]  a  complicated  finite  element  model  was  validated  and  the  natural 
frequencies  were  found  for  various  boundary  conditions.  Those  natural  frequencies  were 

used  as  reference  for  the  selection  of  the  appropriate  thickness  and  modulus  of  elasticity 
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using  the  same  boundary  conditions.  The  natural  frequencies  of  the  ramp  model  were 
calculated  by  solving  the  eigenvalue-eigenvector  problem 


([K]-co^[M]){d}  =  (),  (70) 

where  O)  is  the  natural  frequency  in  radians  per  second  and  jt/}  is  the  corresponding 

mode  shape.  Different  ramp  thickness  and  Young’s  moduli  were  tested,  and  the 
corresponding  natural  frequencies  were  found.  The  natural  frequencies  obtained  from  the 
solution  of  the  above  problem  can  be  seen  in  the  Table  1. 


From  Ref,  [6] 

E  =  39*10^  psi, 

E  =  29*10®  psi, 

E  =  40*10®  psi, 

E  =  49*10®  psi. 

case  5 

t  =  10  in 

t  =  20  in 

t  =  20  in 

t  =  25  in 

2.1 

0.37 

0.639 

0.75 

0.65 

2.62 

1.22 

2.1 

2.47 

2.64 

3.45 

1.34 

2.28 

2.68 

3.45 

6.65 

2.6 

4.46 

5.24 

6.04 

Table  1 .  Natural  Frequencies  of  the  Ramp  Model. 


From  the  above  results  it  is  obvious  that  we  cannot  model  exactly  a  complicated 
structure  like  the  CAPE-D  ramp  as  a  plate,  with  uniform  properties.  Detailed  finite 
element  models  are  more  appropriate  for  that  purpose.  Since  the  purpose  of  this  work  is 
to  establish  the  interaction  between  the  ramp  and  the  RRDF  and  not  to  find  the  exact 
response  and  stress  values  of  the  ramp,  the  plate  model  can  approximate  with  acceptable 
accuracy  the  physical  ramp.  The  only  exception  to  this  is  the  first  natural  frequency, 
which  corresponds  to  longitudinal  bending.  A  closer  look  at  the  physical  structure  of  the 
ramp  reveals  that  the  ramp  is  highly  reinforced  longitudinally  by  stiffeners  that  run 
throughout  the  whole  length  of  the  ramp.  This  explains  the  difference  in  the  first  natural 
frequency  in  bending.  Based  on  the  above  observations,  we  select  the  last  case  presented 
in  Table  1  as  the  baseline  ramp  for  the  current  work. 
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IV.  COUPLING  BETWEEN  SHIP  RAMP  AND  BARGE  MODELS 


A,  INTRODUCTION 

In  section  II  the  mathematical  models  of  the  ship  and  the  barge  in  regular  waves 
were  created.  In  section  III  the  finite  element  of  the  ramp  was  created  and  validated. 
From  those  models  the  mass,  stiffness,  and  damping  matrices  of  each  model  can  be 
found.  Since  damping  was  not  considered  for  the  ramp,  its  damping  matrix  is  zero. 
Previous  studies  (see  for  example  Ref  [5]),  have  suggested  that  active  or  passive  motion 
compensation  can  reduce  the  motions  and  the  stresses  of  the  ramp.  One  such  coupling 
condition  for  motion  and  stress  reduction  is  the  use  of  a  spring-damper  combination  at  the 
two  ends  of  the  ramp,  between  the  barge  and  the  ramp.  This  coupling  method  is  used  in 
the  current  work.  The  coupling  methodology  is,  however,  general  enough  so  that  other 
coupling  conditions  can  be  also  used.  A  sketch  of  the  system  can  be  seen  in  Figure  7. 

B,  COUPLING  BETWEEN  SHIP,  RAMP,  BARGE 

The  equations  of  motion  for  the  ship,  the  ramp  and  the  barge  in  the  time  domain 
are 

[M,  +  4]{4)  +  l5,]fe)  +  [C,]{l7,)={^'„)  (VI) 

[M,]{rf,}  +  [C,]{d,}  +  [i:,]K)  =  {F,)  (72) 

[M,  +  4]{4)  +  [Vi»]{4}  +  [Q]{'74={^',.»)  (73) 

where  the  subscripts  s,  b,  and  p  correspond  to  the  ship,  the  barge,  and  the  ramp 
respectively.  The  matrices  in  the  ship  and  barge  equations  were  defined  in  section  II  and 
the  matrices  in  the  ramp  equation  were  obtained  from  the  finite  element  model  as 
described  in  section  III.  The  matrices  and  represent  the  wave  forces  acting  on 
the  barge  and  on  the  ship.  Note  that  the  damping  matrix  for  the  ramp  is  zero.  Equations 
(71)  and  (73)  have  been  solved  as  a  system  of  twelve  equations  with  twelve  unknowns  in 
section  II. 

Since  the  ship  is  very  heavy  compared  to  either  the  ramp  or  the  barge,  its  inertia 
forces  are  orders  of  magnitude  greater  than  the  barge’s,  it  is  reasonable  to  assume  that  its 
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motions  are  not  going  to  be  affected  by  the  dynamic  forces  transmitted  through  the  ramp. 
Therefore,  the  motions  of  the  ship  can  be  considered  as  uncoupled  from  the  ramp  and  the 
barge  and  they  are  going  to  be  given  by  the  solution  obtained  in  section  II. 

1.  Spring-Damper  Coupling  Condition 

The  main  coupling  condition  considered  in  this  work  is  the  spring-damper  case. 
Between  the  ramp  and  the  barge  we  assume  the  existence  of  two  sets  of  spring-damper 
combinations  attached  to  the  end  points  of  the  ramp,  as  can  be  seen  in  Figure  7  in 
Appendix  B.  From  the  barge’s  equations  of  motion,  the  heave,  roll  and  pitch  equations 
are  going  to  be  affected  from  the  coupling  like  follow: 

Heave  equation: 


^316  ^32*  ^^7,b  ^33*  ^ 34b  ^34b  ^ 35b  ^35b  ^ 36b  ^36b  )  1 


^b 

Vib 

V3b 

V4b 

Vbb 

iVbb 


+ 


(^316  ^32b  ^33b  ^34b  ^35b  ^36b) 


Vib 

nxb 

Vib 

Vib 

V3b 

"“'“(Ql*  ^31b  ^33b  ^34b  ^35b  Qei  ) " 

V3b 

V4b 

V4b 

Vsb 

Vsb 

^6b^ 

^6b. 

{^3wb)  {‘^A  ^a)^A  ^b)^B  {‘^A  ^a)^A  {‘^B 


the  right  hand  side  becomes 

LHS  =  (74) 

where  ^  was  defined  in  Equation  (9)  and  k^,kg,c^,Cg  are  the  spring  and  damping 
coefficients  at  points  A  and  B  ,  respectively.  By  substituting  the  value  of  ^  in  Equation 
(74)  we  get 
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(75) 


LHS  —  {k^  +  )  VAb  {^AyA  ■*■  ^syB  )  ^5b  (  ^A^A  ) 

~  i^A  “*“  )  ~  ^Ab  (^AyA  “*“  )  ~  ^5b  (~^A^A  ~  ^B^B  )  “*“ 

+  w^k^  +  Wj,kg  +  +  w^Cj, 

Equation  (75)  can  be  written  as 

{^ib  ^3*  }  {%  }  {^ihp  }  {fib  }  {^ibp  }  {Vb  }  ~  '^A^A  ~  ^B^B  ~  ^A^A  ~  ^B^B  ~  {^3wb  )  (^6) 

where 

{^3Ap}  “  {^31  ^32  ^33  ^A  ^B  ^34  “*■  ^35  ~  ^A^A  ~  ^B^B  ^ib}  C^'^) 

and 

{^3Ap}  “{^31  ^32  ^33  '^^A'^^B  ^34  ■*■  ■*"  ^35  ~  ^A^A  ~^B^B  ^3b\  (^^) 


Roll  Equation: 


{^A\b  “*“  ^41*  -^^42*  “*“  ^A2b  ^A3b  “*“  ^A2b  ^AAb  “*“  ^446  ^Abb  “*“  ^Abh  ^ Abb  “*“  ^Abb  )  j 


V,b 

bib 

bib 

bAb 

bbb 

bbb 


+ 


( -^41*  ^Alb  ^A3b  ^AAb  ^Abb  ^Abb  ) 


'b,b' 

bxb 

bib 

bib 

bib 

"  +  (^'41^  C^2b  ^Aib  ^AAb  ^Abb  ^Abbf 

bib 

bAb 

bAb 

bbb 

bbb 

bbb  ^ 

bbb. 

=  {FA^b)-{^A-WA)kAyA-{^B-WB)  ^ByB  “  (L  “  %  ) 


the  right  hand  side  becomes 
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LHS  —  )  ^A^Ay A  ‘^B^sys  ^A^AyA  ^B^By B  ■*■ 

+  ^A^AyA  +  ^B^SyB  +  ^A^AyA  +  ^B^Bys 


(79) 


After  massaging  equation  (79)  beeomes 

{M 4,  +  +  {Babp  }  {%  }  +  {C^bp  }  {%  }  -  -  Wgkpyj,  -  w^c^y^  -  WpC^y^ 


where 


(-^4wA  ) 


Kb]  = 


(80) 


{C'41  C42  C43+A:^_y^  k^yg  C^^+k^y^  k^y^  C45  k^x^y  kgX^y  g 


(81) 


and 


Kb]  = 

{-^41  B^2  ^Ay  a  ~  ^By  b  B^i^  +  c^y  ^  —  Cgy  g  B^^—c^x^y^—CgXgyg  B^^ 


(82) 


Pitch  Equation; 


^^5\b~^^\b  ^52b~^^52b  ^52h~^^52h  ^  54b  ^54h  ^  55h  ^55b  ^  56b  ^56b')\ 


Ixb 

Vib 

%b 

V4h 

V,b 

Vih 


(^51Z»  ^52b  ^52h  ^54b  ^55b  ^56b)\ 


(^5wi  )  +  (#4  -  +  (#5  -  )  kgXg  +[^^-W^)c^X^+(^g-Wg)  CgXg 

Easily  ean  be  shown  that  equation  the  above  equation  beeomes 


Vib' 

Vib 

nib 

%b 

V4b 

"■*■(^16  6^ 52b  6^ 52b  ^4*  ^56  ^6*  ) ' 

V2b 

V4b 

Vib 

V5b 

%t\ 

Vbb. 
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{Msh+Ab]{'nb]  +  [Bsbp}{nb}+[Csbp]{nt]+WAkA^A+WBkBXB+W^C^X^+W^C^X^ 

=  [Fs^b ) 


(83) 


where 


K.}  = 

{^51  Q2  —  kgXg  —  k^x^y  ^  —  kgX^y  g  C^^+k^x^+kgXg  Cjgj 


(84) 


and 


{Fsbp  } 

{^51  ^52  —  c^x^  —  CgXg  —  c^x^y  ^  —  CgXgy  g  +  c^x^  +  CgXg  B 


(85) 


For  the  ramp  the  equation  of  motion  in  the  time  domain  beeomes 

[_^  p\[^p]^[Fp~\[dp]  +  [Fj\[dp]  =  [Fp]  +  {^^-w^)k^+{^g-Wg)kg  + 

+  (L-Wa)ca  +  (L-Wb)cb  (86) 

Substituting  ^  and  rearranging,  Equation  (86)  beeomes 

[Mj\[dp]  +  [cJ\[d^]  +  [Kj\[d^]-%{k^  +  kg)-ri,{k^y^+kgyg)- 

-Vb  {-^A^A  -  ^B^B  )-%{Ca+Cb)-  %  {cAyA  +  "  Vs  (-^A^A  “  )  + 

+  W/^  +  Wgkg  +  +  WgCg  =  [F ^  }  (87) 


or 

[_^pWp]^\.^pb\{^p]^  V^pb  ~\{dp]-'nAkA+kB)-VA  {^AyA + ^ByB )  - 

-%  {-Xa^a  -  ^BkB  )-vAca+Cb)-  %  (c^yj  +  Cgyg )  -  %  (-x^c^  -  XgCg )  =  {e;}  (88) 

where 
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C„ 


c 


363  1 


Q33  333 


Q 


363 


■■  C  +c 

'-'363  363  ^  '-'A 


and 


(89) 


■  ^11 

v 

363 

•  •  •  fC  -\-  k  •  •  • 

^^333  333 

^363  I 

...  V 

^^363  363  ^ 

2,  Roll  Mitigation  Coupling  Condition 

The  second  coupling  condition  considered  is  a  roll  mitigation  coupling  condition. 
For  demonstration  purposes,  we  assume  that  one  point  of  the  barge  is  connected  at  the 
center  of  the  ramp  with  the  use  of  a  special  universal  joint.  This  allows  heave  transfer 
from  the  barge  to  the  ramp  but  leaves  the  two  other  degrees  of  freedom  (roll  and  pitch) 
free.  The  universal  joint  will  be  modeled  as  a  spring  that  connects  one  point  at  the  end  of 
the  barge  with  three  points  in  the  middle  of  the  ramp.  The  three  points  on  the  ramp  were 
considered  to  represent  an  area  in  the  physical  model.  The  stiffness  of  the  springs  was 
chosen  such  that  the  vertical  motion  of  the  ramp  follows  the  vertical  motion  of  the  barge 
without  producing  singularities  in  the  solution. 


By  constructing  a  free  body  diagram  for  the  ramp  and  the  barge  the  new  equations 
of  motion  can  be  found.  The  coupled  equations  of  motion  of  the  barge  for  the  heave,  roll, 
pitch  motions  become 


+  +  ={^'«}-£(#c  -»’c,  h  (91) 

+  {S«}{'7»}  +{C«}{)7,}  =  {F,.*}  + 1  (4  -  We,  (92) 
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{M,,  +  +{^3  Ji'?.} = {^;.*} + 2;(4  -  «'c.  ]Kd„  (93) 

i-1 

where  A:^  is  the  stiffness  of  the  springs,  and  C,C,  represent  the  conneetion  points  on  the 
barge  and  the  ramp  respectively.  By  substituting  the  value  of  ^  in  Equations  (91)  -  (93) 
and  rearranging  we  get 


[M,  +A]{%]  +  [Bt]{%]  +  [Ckp][%]  = 


where 


F, 


\wb 


2wb 

3 


^3wb  ^C, 

i-1 

3 

^4wb  ^cd 


i=\ 

3 


d^Swb  ^  ^  dx, 


i=\ 


F. 


6wb 


(94) 


K] 


c  c 

'-'11  '-'12 

c  c 

'-'21  '-'24 


c  c 

''"31  '-'32 


c 


c 


23 


^^3  3  ^  ^ 


i=\ 

3 


c, 


14 


c 


24 


^4  ^  dyq  k(- 

/-I 

Q,  Q  Qj+Z4k,*c  C44+E‘'L*c 


(=1 


c, 


15 


c 


25 


(=1 


t=l 


/=! 


/=1 


/=1 


c  c 

^61  ^62 


c 


63 


c 


64 


c 


65 


For  the  ramp  the  equation  of  motion  in  the  time  domain  becomes 

(=1 

By  substituting  ^  and  rearranging,  Equation  (96)  becomes 


C, 


16 


c 


26 


D 

^35  ~  ^  ^xC,.  6'36 


6*45  ^  ^xc,  Byq 


6*51  ^'52  6'53  C54  ^  ^xc,  ^yc,.  6*55 -I- ^  C^g 


C 
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(95) 


(96) 
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or 


_^p]{'^p}+[<^JR}+MR}  +3%^c^c  +  &q^c  =R}  (97) 


[Mj\{dp]  +  [C,]{d^]  +  [^t:^,~\[d^]-3%kc-3ri,yckc+3ri^Xckc={F^}  (98) 

where 


■  ^11 

-^1  363 

•••  K..+kc  •••  : 

^363  1 

^363  363 

and  i  =  180,183,186. 

C.  COUPLED  SOLUTION  IN  FREQUENCY  DOMAIN 

The  solution  of  the  ship-ramp-barge  system  in  the  frequency  domain  is  required  in 
order  to  obtain  the  new  RAO  of  the  barge  motions  and  the  RAO  of  the  ramp  motion.  In 
order  to  obtain  the  solution  of  the  complete  system,  the  motion  of  the  ship  has  to  be  taken 
into  account.  It  is  assumed  that  the  end  of  the  ramp  attached  to  the  ship  is  pinned  on  the 
ship’s  stern.  Therefore,  all  the  points  along  the  ship-end  of  the  ramp  are  going  to  follow 
the  ship’s  vertical  motion.  This  induces  the  following  boundary  condition  for  all  points 
across  the  ship-end  of  the  ramp 

«',=£,  (100) 

and 

(».),=-%,  (101) 

where  w,  is  the  displacement  along  the  z  axis  of  a  point  i  on  the  ramp,  {0^).  is  the 
rotation  along  the  x  axis  of  a  point  i  on  the  ramp,  is  the  displacement  along  the  z  axis 
of  a  point  i  on  the  ship,  and  is  the  roll  angle  of  the  ship. 
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Quantities  ^,,74  are  known  from  the  hydrodynamies  of  the  ship.  Therefore, 

Equations  (100)  and  (101)  ean  be  implemented  as  boundary  eonditions  in  the  finite 
element  model  of  the  ramp.  This  means  that  the  stiffness  and  the  foree  matriees  of  the 
ramp,  as  shown  in  Equations  (90)  and  (99)  have  already  taken  into  aeeount  the  ship’s 
motions. 

It  should  be  pointed  out  that  the  equations  of  motion  of  the  barge  as  shown  in  (73) 
do  not  refieet  the  proximity  effeet  of  the  ship’s  presenee.  To  take  into  aeeount  the  ship’s 
presenee,  additional  forees  have  to  be  added.  In  the  solution  presented  in  section  II  those 
forces  were  taken  into  account  since  the  full  12x12  system  was  solved,  while  Equations 
(73)  are  only  a  6x6  system.  Eorces  due  to  proximity  effects  are  added  in  Equations  (76) 
and  (94)  as  additional  external  excitation  vectors  {fyj}  that  are  proportional  to  the  ship’s 
accelerations  and  velocities. 

Eurthermore,  in  the  scenario  considered  in  this  work,  the  load  of  two  tanks  was 
placed  in  the  middle  of  the  ramp.  This  was  done  in  order  to  model  the  envisioned  worst- 
case  scenario  of  one  vehicle  stopped  because  of  mechanical  failure  and  another  coming  to 
provide  assistance.  This  adds  an  additional  force  vector  in  Equations  (88)  and  (98). 
Therefore,  the  equations  of  the  ramp  and  the  barge  can  be  coupled  together  in  a  single 
system  of  equations 


where 


K]  [0] 

[0]  [m,  +  4] 


[c,d  [c.l 
[Q]  K] 


(103) 


(104) 


(105) 
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(106) 


[^2]  [c.,] 


(107) 


and  Cj ,  Cj ,  the  coupling  matrices,  created  by  the  additional  left  hand  term 

quantities  in  the  Equations  (76),  (88),  (94)  and  (98). 

By  using  the  same  linear  theory  assumption  as  in  section  II  the  vector  of 
displacements  can  be  written  as 

(109) 


where  again  0)^  is  the  frequency  of  the  incoming  wave  and  are  the  complex 
amplitudes  of  the  displacements. 

The  solution  of  Equation  (102)  in  the  frequency  domain  then  becomes 


[Fgu,] 


(110) 


where  are  the  complex  amplitudes  of  the  forces. 

The  above  solution  represents  the  Response  Amplitude  Operator  (RAO)  of  the 
ramp  and  the  barge  in  the  coupled  system. 

D.  STRESS  CALCULATION  FOR  THE  RAMP 

The  solution  of  Equation  (110)  can  be  used  to  obtain  the  stress  for  the  ramp.  Erom 
Hooke’s  Eaw  it  is  known  that 


a  =  E£ 

(111) 

II 

(112) 

Substituting  Equations  (40),  (41)  in  (1 1 1),  (1 12)  and  rearranging 


40 


(113) 


~E 

0 

o: 

>  = 

0 

E 

0 

< 

1  1 

0 

0 

G\ 

/ - 

The  same  proeess  with  Equations  (42),  (43)  yields 


(114) 


The  strains  can  be  found  by  substituting  the  solution  of  Equation  (110)  in 
Equations  (61)  and  (62).  Therefore,  the  Von-Misses  stress  for  the  elements  of  the  ramp 
can  be  found  from 


'VM 


2  2  2 

+  r  +r 

xy  yz  xz 


(115) 


where  in  our  case  <7^=0.  Unlike  the  ramp  validation  case,  the  true  value  of  Young’s 
modulus  of  elasticity  is  used  to  calculate  stresses. 

E.  RESULTS 

Numerical  results  are  presented  in  Appendix  B.  The  general  parameters  used  in 
the  numerical  solution  are  presented  in  Table  2.  The  spring-damper  boundary  conditions 
were  used  throughout  the  results,  see  Eigure  7  of  Appendix  B  for  a  schematic.  The 
nominal  values  for  the  spring  and  damper  constants  are  also  shown  in  Appendix  B. 
Unless  otherwise  stated,  these  values  will  be  used  for  the  results. 

The  effect  of  the  flexural  rigidity  of  the  ramp  and  the  assumed  boundary 
conditions  on  RRDE  response  is  shown  in  Eigure2  B1  through  B24.  It  can  be  seen  that,  in 
general,  the  effects  of  the  ramp  does  not  alter  the  fundamental  frequency  response  shape 
of  the  RRDE.  This  can  be  attributed  to  the  high  hydrostatic  stiffness  of  the  large 
waterplane  area,  shallow  bodies,  that  make  up  the  RRDE.  Although  the  general  shape  of 
the  response  amplitude  operator  does  not  change,  the  relative  values  differ  so  that  the 
response  in  a  seaway  will  be  different.  This  means  that  the  effect  of  the  ramp  and  the 
spring/damper  constants  on  the  RRDE  is  expected  to  be  measurable  and  it  could  be 
quantified  in  random  seas.  This  will  be  explored  further  in  the  next  chapter. 
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A  comprehensive  set  of  plots  showing  the  response  of  the  RRDF  as  modified  by 
the  ramp,  for  all  heading  angles  is  presented  in  Figures  B25  through  B30.  As  expected, 
wave  direction  has  a  very  signifieant  effect  on  RRDF  response.  Therefore,  results  in 
random  seas  are  expected  to  vary  greatly  depending  on  sea  direetionality. 

Figures  B31  through  B38  present  results  on  the  relative  motions  between  at  the 
two  touchdown  points  between  the  ramp  and  the  RRDF.  Point  A  is  the  starboard 
connection  point,  and  point  B  is  the  port  connection  point,  see  Figure  7.  Figures  B39  and 
B40  summarize  the  relative  motion  results  for  all  wave  headings.  The  relative  motion 
results  are  used  in  figures  B41  and  B42  to  calculate  and  present  the  ramp  twist  angle 
RAO  as  a  function  of  frequency  and  wave  direetion.  The  ramp  twist  angle  is  of  extreme 
importance  in  operations,  since  it  is  directly  related  to  maximum  stress,  whieh  is  in  most 
cases  the  limiting  factor  in  performance.  Figures  B43  through  B46  show  3-dimensional 
views  of  the  ramp  displacement  field.  These  are  presented  seleeted  frequencies  and  wave 
directions  and  are  to  be  used  for  general  visualization  purposes;  the  vertieal  scale  is  mueh 
smaller  than  the  horizontal  seale  and  do  not  necessarily  refleet  the  aetual  ramp 
deformation  in  a  forced  motion  setting. 

Finally,  the  remaining  figures  in  Appendix  B  show  representative  RAO’s  in  terms 
of  ramp  stress.  The  maximum  stress  in  the  ramp  in  terms  of  frequency  and  wave  direetion 
is  shown  in  Figures  B47  and  B48.  Peaks  in  the  stress  eorrespond  to  either  maximum 
relative  motion  or  ramp  twist.  When  coupled  with  an  appropriate  seaway  model,  these 
will  result  in  maximum  expected  values  in  a  given  seaway.  Representative  stress  fields  in 
the  ramp  for  select  wave  periods  and  headings  are  shown  in  Figures  B49  through  B52. 
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V.  RESPONSE  IN  A  SEAWAY 


A,  INTRODUCTION 

The  ultimate  goal  of  ship  motion  studies  is  to  be  able  to  prediet  how  the  ship  will 
behave  in  realistie  irregular  seas.  As  already  mentioned  in  seetion  II  the  ealculation  of  a 
vessel’s  motions  in  regular  waves  is  only  the  first  step  towards  the  prediction  of  the 
motions  in  a  seaway.  According  to  the  method  proposed  in  1953  by  St.  Denis  and 
Pierson,  using  spectral  methods  developed  in  other  fields,  it  is  possible  to  relate  the 
spectral  density  of  ship  responses  to  the  input  ocean  wave  spectrum.  In  Ref  [2]  and  [8] 
the  interested  reader  can  find  a  very  detailed  discussion  on  seaway  modeling.  Here  only  a 
short  description  and  some  useful  formulas  will  be  presented. 

B,  DESCRIPTION  OF  THE  SEAWAY 

The  wave  patterns  in  an  open  sea  are  ever  changing  with  time  and  space,  in  a 
manner  that  appears  to  defy  analysis  be  it  linear  or  second  order  Stokes  waves.  Ambient 
waves  on  the  surface  of  the  sea  are  dispersive  as  well  as  random.  Random  refers  to  the 
character  of  the  wave  height  distribution.  The  generating  mechanism  is,  predominantly, 
the  effect  upon  the  water  surface  of  wind  in  the  atmosphere.  The  wind  is  itself  random, 
especially  when  viewed  from  the  standpoint  of  the  turbulent  fluctuations  and  eddies 
which  are  important  in  generating  waves.  The  randomness  of  sea  waves  is  subsequently 
enhanced  by  their  propagation  over  large  distances  in  space  and  time  and  their  exposure 
to  the  random  nonuniformities  of  the  water  and  air. 

According  to  the  principle  of  superposition,  one  could  attempt  a  description  of  the 
seaway  deterministically  accounting  for  each  one  of  the  individual  wave  components. 
Leaving  aside  the  issue  of  whether  such  an  approach  is  possible  or  not,  it  is  clearly  a 
highly  impractical  task.  Fortunately,  it  is  not  required  either.  It  would  suffice  to  represent 
sea  waves  in  a  probabilistic  manner. 

The  statistical  properties  of  a  random  process  can  be  defined  with  two  possible 
ways.  It  can  be  considered  with  statistical  properties  taken  "across  the  ensemble"  at  fixed 
values  of  iimQt  =  t^,t  =  1^,  etc.,  or  may  be  considered  with  properties  of  the  random 
process  taken  "along  the  ensemble"  where  etc.,  are  assumed  to  vary.  Description  of 
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a  random  process  so  general  would  require  an  enormous  amount  of  information. 
Fortunately,  for  certain  random  processes  such  as  sea  waves,  it  is  possible  to  assume  that 
the  proeess  is  of  a  special  form:  "stationary",  "homogeneous",  and  "ergodic".  For  sea 
waves  it  should  be  adequate  to  describe  the  wave  environment  over  a  period  of  a  few 
hours  (or  before  the  next  weather  report  comes  in),  and  to  assume  that  the  wave  motion  is 
stationary  (its  statistics  remain  the  same  over  time)  during  this  interval  of  time.  Likewise, 
there  is  no  interest  in  describing  the  wave  environment  simultaneously  throughout  the 
world,  a  small  area  of  operations  is  all  is  needed,  and  it  can  be  assumed  that  the  wave 
motion  is  homogeneous  (its  statistics  remain  the  same  in  that  area)  in  space  over  the  area 
in  question.  These  statements  have  meaning  only  in  a  statistical  sense,  since  for  a  random 
wave  system  it  would  be  funny  to  suggest  that  the  precise  wave  motion  is  the  same  at 
different  points  in  space  or  time.  Eventually,  it  is  reasonable  to  assume  that  the  statistical 
properties  of  the  waves  measured  over  time,  should  be  "typical"  of  the  random  process. 
This  means  that  they  should  be  the  same  even  if  we  were  able  to  sample  all  possible 
realizations  of  the  wave  motion  at  a  fixed  time.  Such  random  processes  are  said  to  be 
ergodic. 

One  of  the  most  significant  parameters  needed  in  order  to  arrive  at  a  statistical 
description  of  the  seaway  is  the  total  mean  energy  of  the  wave  system  per  unit  area  of  the 
free  surface  (E ).  It  can  be  shown  from  hydrodynamics  that  this  is  equal  to: 

S{co,0)d6dco  (116) 

where  6  is  the  wave  direction.  The  function  S[co,6)  is  called  the  spectral  density  and 

more  rigorously  can  be  defined  as  the  Fourier  transform  of  the  correlation  function  for  the 
free  surfaee  elevation.  It  is  customary  to  ignore  the  factor  pg  and  to  refer  to  the  function 

S[co,6)  as  the  spectral  energy  density  or  simply  the  energy  speetrum.  More  speeifically, 

this  is  a  direetional  energy  spectrum.  It  can  be  integrated  over  all  wave  directions  to  give 
the  frequency  spectrum: 

=  S  {(0,0)  dO 

Using  Equations  (116)  and  (1 17)  we  can  see  that: 
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(117) 


(118) 


£■  =  £  S[co)dco 

i.e.,  the  area  under  the  speetrum  S{co)  is  equal  to  (within  the  constant  pg)  the  mean 

energy  stored  in  that  particular  wave  system.  If  one  attempts  to  find  the  sea  wave 
spectrum  from  measurements  of  the  free  surface  elevation  at  a  single  point  in  space,  for 
instance  by  recording  the  heave  motion  of  a  buoy,  the  directional  characteristic  of  the 
waves  will  be  lost.  Only  the  frequency  spectrum  (117)  can  be  determined  from  such  a 
restricted  set  of  data.  A  limited  amount  of  directional  information  follows  if  one  measures 
the  slope  of  the  free  surface,  for  example  by  measuring  the  angular  response  of  the  buoy 
as  well  as  its  heave.  A  complete  description  of  the  directional  wave  spectrum  requires  an 
extensive  array  of  measurements  at  several  adjacent  points  in  space,  and  there  are 
practical  difficulties  associated  with  this  task. 

As  a  simpler  alternative,  one  can  assume  that  the  waves  are  unidirectional,  with 
the  energy  spectrum  proportional  to  a  delta  function  in  0 .  Wave  spectra  of  this  form  are 
called  long  crested,  since  the  fluid  motion  is  two-dimensional  and  the  wave  crests  are 
parallel,  and  the  frequency  spectrum  (1 17)  is  sufficient  to  describe  the  wave  environment. 

If  the  waves  are  generated  by  a  single  storm,  far  removed  from  the  point  of 
observation,  it  might  be  presumed  that  these  waves  would  come  from  the  direction  of  the 
storm  in  a  long  crested  manner.  The  limitations  of  this  assumption  are  obvious  to  anyone 
who  has  observed  the  sea  surface.  While  a  preferred  direction  may  exist,  especially  for 
long  swell  that  has  traveled  large  distances,  even  these  long  waves  will  be  distributed  in 
their  direction,  and  for  short  steep  waves  the  directional  variation  is  particularly 
significant.  Since  the  superposition  of  such  waves  from  a  range  of  different  directions 
appears  in  space  as  a  variation  of  the  free  surface  elevation  in  all  directions,  these  waves 
are  known  as  short  crested  waves. 

Usually  in  the  fields  of  ocean  engineering  and  naval  architecture  it  is  customary  to 
assume  that  the  waves  are  long  crested.  With  such  a  simplification  it  is  possible  to  use 
existing  information  for  the  frequency  spectrum  (117),  which  is  based  on  a  combination 
of  theory  and  full-scale  observations.  Sea  wave  spectra  depend  on  the  velocity  of  the 
wind  as  well  as  its  duration  in  time  and  the  distance  over  which  the  wind  is  acting  on  the 
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free  surface.  This  distance  is  known  as  the  fetch.  Wave  spectra  that  have  reached  a  steady 
state  of  equilibrium,  independent  of  the  duration  and  fetch  are  known  as  fully  developed. 
There  are  many  proposals  of  the  analytical  expression  of  the  frequency  spectrum.  Some 
of  them  are  generic  and  some  of  them  are  based  on  observations  and  data  collected  in 
specific  geographical  regions.  A  semi-empirical  generic  expression  for  the  frequency 
spectrum  of  fully  developed  seas  is: 


5(6;) 


(O 


exp 


f 

Uo) 


(119) 


Here  a  and  (5  are  nondimensional  parameters  defining  the  spectrum,  and  U  is  the  wind 
velocity  at  a  standard  height  of  19.5  meters  above  the  free  surface.  This  two  parameters 
spectrum  is  sufficiently  general  to  fit  quite  a  few  observations  and  is  consistent  with 
theoretical  predictions  of  the  high  frequency  limit.  The  most  common  values  for  these 
parameters  are: 

a  =  8.1x10"'  and  J3  =  0.74  (120) 

and  with  these  values  it  is  known  as  the  Pierson  -  Moscowitz  spectrum. 

C.  SHIP/RAMP/BARGE  RESPONSE  SPECTRAL  ANALYSIS 

The  relationship  between  motions  in  regular  and  irregular  waves  must  be 
established  in  order  to  use  either  theoretical  equations  or  experimental  studies  to  lead  to 
practical  predictions  of  system  behavior  in  a  realistic  seaway.  The  required  relationship  is 
expressed  by  the  principle  of  linear  superposition,  which  says  that  the  response  of  a  ship 
to  an  irregular  sea  can  be  represented  by  a  linear  summation  of  its  responses  to  the 
component  regular  waves  that  comprise  the  irregular  sea.  To  make  use  of  this  principle,  a 
ship  response  to  regular  waves,  whether  it  has  been  determined  by  calculation  or  by 
model  measurement,  is  expressed  in  a  Response  Amplitude  Operator  (RAO)  form  as 
measure  of  the  response  to  a  regular  wave  of  unit  amplitude. 

Like  the  sea  waves  themselves,  a  ship  response  is  a  random  variable.  The 
statistics  of  a  floating  body  response  are  identical  to  the  wave  statistics,  except  that  the 
wave  energy  spectrum  S  is  multiplied  by  the  square  of  the  RAO  (this  is  a  property  of 
linear  systems).  Thus,  if  the  subscript  R  represents  any  body  response,  we  have: 
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Sj,{cD)  =  \z^{cDf  S{0))  (121) 

where  Zj^[co)  is  the  RAO  of  the  response  R,  and  the  speetrum  of  the  seaway. 

Equation  (121)  can  then  be  utilized  to  obtain  the  spectrum  of  the  response  R  ,  and  after 
that  specific  formulas  can  be  used  to  provide  the  statistics  of  this  response  in  a  particular 
sea  state. 

To  a  large  extent,  Equation  (121)  provides  the  justification  for  studying  regular 
wave  responses.  The  transfer  function  [co]  is  valid  not  only  in  regular  waves,  where  it 

has  been  derived,  but  also  in  a  superposition  of  regular  waves,  and  ultimately  in  a 
spectrum  of  random  waves.  Generally  speaking,  a  vessel  with  favorable  response 
characteristics  in  regular  waves  will  be  good  in  irregular  waves,  and  vice  versa.  This 
statement  is  oversimplified,  however,  and  the  relative  shape  of  the  energy  spectrum  and 
the  transfer  function  is  very  crucial.  Eor  example,  a  large  resonant  response  of  the  body 
will  be  of  importance  only  if  the  resonant  frequency  is  located  close  to  the  peak  of  the 
wave  energy  spectrum,  and  vice  versa. 

Eventually,  the  statistical  predictions  of  the  amplitudes  of  ship  motions  may  be 
determined  from  the  variance,  or  spectral  area,  of  a  ship  response  spectrum  given  by  the 
equation: 

(122) 

J  0 

In  order  to  present  the  random  wave  results  in  a  more  compact  manner,  the 
significant  and  the  Root  Mean  Square  (RMS)  values  of  the  responses  will  employed.  The 
significant  value  of  a  random  process  is  defined  as  the  average  of  the  1/3  of  the  highest 
values  of  all  responses  and  is  frequently  used  in  design  to  characterize  the  severity  of  a 
particular  response.  The  RMS  value  can  be  easily  computed  from  the  spectrum  of  the 
response  and  is  in  the  square  root  of  the  area  under  the  spectrum  of  the  response,  which  is 
the  integral  of  the  spectrum  function  over  the  frequency  range.  The  significant  value 
(single  amplitude)  is  equal  to  two  times  the  RMS  value.  To  get  the  “8  hour  maximum 
amplitude”,  the  RMS  vale  must  be  multiplied  by  four.  This  is  also  equal  to  the  significant 
double  amplitude  or  “height”  (peak  to  trough  deviation)  of  the  response. 
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D,  RESULTS 

Typical  results  in  random  seas  are  presented  in  Appendix  C.  For  demonstration 
purposes,  we  use  the  standard  Pierson-Moskowitz  speetrum  at  sea  state  3,  shown  in 
Figure  Cl.  This  corresponds  to  fully  developed,  long  erested  unidireetional  seas  of 
approximately  4  ft  significant  wave  height,  significant  range  of  periods  between  2  and  7 
seconds,  an  average  period  of  4  seconds,  a  wind  speed  of  15  knots,  and  waves  of 
approximately  50  ft  wave  length.  We  also  use  three  different  settings  of  the  spring  and 
damper  coeffieients,  as  follows: 

•  (A=l,  c=l)  eorresponds  to  a  soft  spring/damper  combination  with  both  values 
equal  to  the  nominal  values  shown  in  Table  2  of  Appendix  B. 

•  {k=2,  c=2)  eorresponds  to  a  hard  spring/damper  combination  with  both  values 
equal  to  10  times  the  nominal  values  for  k  and  c. 

•  (A=3,  c=3)  corresponds  to  a  very  hard  spring/damper  eombination  with  both 
values  100  times  larger  than  the  nominal  values  shown  in  Table  2  of 
Appendix  B. 

Sample  plots  of  response  speetra  are  shown  in  Figures  C2  through  C13,  in  terms 
of  the  spectrum  of  the  relative  motion  at  conneetion  point  A.  The  various  peaks  in  the 
spectrum  correspond  to  combinations  of  peaks  of  the  eorresponding  frequeney  response 
of  the  relative  motion  and  the  peak  of  the  wave  spectrum. 

The  RMS  value  of  the  relative  motion  at  eonneetion  point  A  is  ealculated  based 
on  the  above  speetra  and  is  shown  in  Figure  C14.  Based  on  these  results  we  ean  observe 
the  following: 

1.  The  RMS  values  are  higher  at  port  beam  and  starboard  beam  seas.  For 
direetional  sea  speetra,  the  eurves  viewed  as  funetions  of  heading  would  be 
smoother. 

2.  Port/starboard  response  is  not  symmetric  as  has  been  diseussed  previously. 

3.  Higher  spring  and/or  damping  coefficients  result  in  less  relative  motions,  as 
expeeted.  This  is  more  pronouneed  in  the  case  of  beam  seas  and  less  in  the 
case  of  following  or  head  seas. 
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Figures  C15  through  C20  present  the  results  of  Figure  C14  in  a  set  of  polar  plot 
contours.  The  radial  coordinate  in  the  polar  plots  is  the  RMS  value  of  the  relative  motion, 
while  the  angular  coordinate  designates  the  wave  heading.  As  can  be  seen,  both  spring 
and  damper  constants  have  a  direct  effect  on  the  contour  plots,  with  the  greatest  effect 
due  to  the  spring  constant.  Figures  C21  through  C27  present  similar  results  in  terms  of 
the  ramp  twist  while  Figures  C28  through  C34  present  results  in  terms  of  the  ramp 
maximum  stress.  Here,  higher  values  of  the  spring  constant  result,  in  general,  in  higher 
levels  of  twist  and  stress.  The  above  sets  of  graphs  can  be  used  to  select  a  compromising 
set  of  values  of  spring  constant  and  damping  ratio,  which  will  minimize  both  relative 
motions  and  stresses  in  a  given  seaway. 

The  last  sets  of  graphs.  Figures  C35  and  C36,  present  for  completeness  the  RMS 
values  of  the  barge  motion  at  points  A  and  B  as  functions  of  heading,  and  parameterized 
in  terms  of  the  spring  and  damper  constants. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A,  CONCLUSIONS 

Based  on  the  methodology  and  the  results  presented  in  this  thesis,  the  following 
eonelusions  ean  be  made; 

1.  The  dynamies  of  the  two  ships  show  signifieant  hydrodynamie  eoupling 
between  vertieal  and  horizontal  plane  motions.  More  speeifically, 
vertical/horizontal  plane  eoupling  for  the  ship  is  small  for  head  and 
following  seas,  but  it  ean  be  signifieant  for  beam  seas. 

2.  Hydrodynamie  eoupling  between  horizontal  and  vertieal  plane  motions  for 
the  RRDF  is  very  signifieant  throughout  the  wave  direction  and  frequency 
range  and  cannot  be  neglected. 

3.  Structural  coupling  of  the  ramp  does  not  alter  the  fundamental  frequency 
response  shape  of  the  RRDF.  However,  although  the  general  shape  of  the 
response  amplitude  operator  does  not  change,  the  relative  values  differ  so 
that  the  response  in  a  seaway  is  different. 

4.  Relative  motions  at  the  ramp/RRDF  connection  points,  as  quantified  by 
their  RMS  values  in  random  seas,  are  higher  for  port  beam  and  starboard 
beam  seas. 

5.  Higher  spring  and/or  damper  coefficients  result  in  smaller  relative  motions 
but  higher  levels  of  twist  and  stresses.  These  results  are  not  always  strictly 
monotonic,  which  suggests  that  optimization  based  on  a  systematic  series 
of  parametric  studies  is  possible. 

6.  Isolator  spring  constant  has,  in  general,  a  more  pronounced  effect  on  the 
RMS  values  of  the  response  than  the  isolator  damping  coefficient. 
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B,  RECOMMENDATIONS 

The  primary  recommendations  for  further  research  are  as  follows: 

1.  A  systematic  parametric  study  of  various  isolator  properties  (spring  and 
damper  constants)  and  ramp  loading  conditions  in  short  crested  seas  [Ref 
8], 

2.  Incorporation  of  more  accurate,  computationally  intensive,  commercial 
finite  element  packages  for  the  ramp,  based  on  the  structural  synthesis 
method  presented  in  [Ref  9]. 
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APPENDIX  A 


This  Appendix  contains  the  following  figures: 

•  Figures  A1  through  A12:  Heave,  Roll,  and  Pitch  Amplitudes  for  the  Ship 
and  the  RRDF  for  Various  Wave  Directions. 

•  Figures  A13  through  A28:  Comparisons  Between  Full  and  Reduced 
Hydrodynamic  Models  for  the  Ship  and  the  RRDF. 

•  Figures  A29  through  A44:  Vertical  Motions  at  the  Ramp/RRDF 
Connection  Points. 

Discussion  of  these  figures  is  presented  in  Chapter  II,  Section  D. 
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Figure  Al.  Plot  of  Heave  Amplitude  per  [ft]  Wave  Height,  for  the  Ship  (CAPE-D)  at 

Various  Wave  Directions 


Figure  A2.  Plot  of  Heave  Amplitude  per  [ft]  Wave  Height,  for  the  Ship  (CAPE-D)  at 

Various  Wave  Directions 
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Amplitude  [ft/ft]  •  Amplitude  [ft/ft] 


Plot  of  Heave  Amplitude  per  [ft]  Wave  Height,  for  the  Barge  (RRDF)  at 
Various  Wave  Direetions 


Figure  A4.  Plot  of  Heave  Amplitude  per  [ft]  Wave  Height,  for  the  Barge  (RRDF)  at 

Various  Wave  Directions 
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Wave  period  T  [sec] 

Figure  A5.  Plot  of  Roll  Amplitude  per  [ft]  Wave  Fleight,  for  the  Ship  (CAPE-D)  at 

Various  Wave  Direetions 
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Figure  A6.  Plot  of  Roll  Amplitude  per  [ft]  Wave  Height,  for  the  Ship  (CAPE-D)  at 

Various  Wave  Direetions 
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Figure  A7.  Plot  of  Roll  Amplitude  per  [ft]  Wave  Height,  for  the  Barge  (RRDF)  at 

Various  Wave  Directions 


Figure  A8.  Plot  of  Roll  Amplitude  per  [ft]  Wave  Height,  for  the  Barge  (RRDF)  at 

Various  Wave  Directions 
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Figure  A9.  Plot  of  Pitch  Amplitude  per  [ft]  Wave  Height,  for  the  Ship  (CAPE-D)  at 

Various  Wave  Directions 


Figure  AlO.  Plot  of  Pitch  Amplitude  per  [ft]  Wave  Height,  for  the  Ship  (CAPE-D)  at 

Various  Wave  Directions 
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Figure  A1 1 .  Plot  of  Pitch  Amplitude  per  [ft]  Wave  Height,  for  the  Barge  (RRDF)  at 

Various  Wave  Directions 


Figure  A12.  Plot  of  Pitch  Amplitude  per  [ft]  Wave  Height,  for  the  Barge  (RRDF)  at 

Various  Wave  Directions 
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Figure  A13.  Comparison  between  Full  12x12  Solution  and  Redueed  Approximation  for 
the  Heave  Amplitude  per  [ft]  Wave  Height  at  Following  Seas  (0  deg). 


Figure  A14.  Comparison  between  Full  12x12  Solution  and  Reduced  Approximation  for 
the  Heave  Phase  Angle  per  [ft]  Wave  Height  at  Following  Seas  (0  deg). 
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Figure  A15.  Comparison  between  Full  12x12  Solution  and  Redueed  Approximation  for 
the  Heave  Amplitude  per  [ft]  Wave  Height  at  Heading  Seas  (180  deg). 


Figure  A16.  Comparison  between  Full  12x12  Solution  and  Reduced  Approximation  for 
the  Heave  Phase  Angle  per  [ft]  Wave  Height  at  Heading  Seas  (180  deg). 
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Figure  A17.  Comparison  between  Full  12x12  Solution  and  Redueed  Approximation  for 
the  Heave  Amplitude  per  [ft]  Wave  Height  at  Starboard  Beam  Seas  (90  deg). 


Figure  A18.  Comparison  between  Full  12x12  Solution  and  Reduced  Approximation  for 
the  Heave  Phase  Angle  per  [ft]  Wave  Height  at  Starboard  Beam  Seas  (90  deg). 
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Wave  period  T  [sec] 

Figure  A19.  Comparison  between  Full  12x12  Solution  and  Redueed  Approximation  for 
the  Fleave  Amplitude  per  [ft]  Wave  Fleight  at  Port  Beam  Seas  (270  deg). 


Figure  A20.  Comparison  between  Full  12x12  Solution  and  Reduced  Approximation  for 
the  Heave  Phase  Angle  per  [ft]  Wave  Height  at  Port  Beam  Seas  (270  deg). 
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Figure  A2 1 .  Comparison  between  Full  12x12  Solution  and  Redueed  Approximation  for 
the  Piteh  Amplitude  per  [ft]  Wave  Fleight  at  Following  Seas  (0  deg). 


Figure  A22.  Comparison  between  Full  12x12  Solution  and  Reduced  Approximation  for 
the  Pitch  Phase  Angle  per  [ft]  Wave  Height  at  Following  Seas  (0  deg). 
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Figure  A23.  Comparison  between  Full  12x12  Solution  and  Redueed  Approximation  for 
the  Pitch  Amplitude  per  [ft]  Wave  Height  at  Head  Seas  (180  deg). 


Figure  A24.  Comparison  between  Full  12x12  Solution  and  Reduced  Approximation  for 
the  Pitch  Phase  Angle  per  [ft]  Wave  Height  at  Head  Seas  (180  deg). 
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Figure  A25.  Comparison  between  Full  12x12  Solution  and  Redueed  Approximation  for 
the  Pitch  Amplitude  per  [ft]  Wave  Fleight  at  Starboard  Beam  Seas  (90  deg). 


Figure  A26.  Comparison  between  Full  12x12  Solution  and  Reduced  Approximation  for 
the  Pitch  Phase  Angle  per  [ft]  Wave  Height  at  Starboard  Beam  Seas  (90  deg). 
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Figure  A27.  Comparison  between  Full  12x12  Solution  and  Redueed  Approximation  for 
the  Piteh  Amplitude  per  [ft]  Wave  Fleight  at  Port  Beam  Seas  (270  deg). 


Figure  A28.  Comparison  between  Full  12x12  Solution  and  Reduced  Approximation  for 
the  Pitch  Phase  Angle  per  [ft]  Wave  Height  at  Port  Beam  Seas  (270  deg). 
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Figure  A29.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  A, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  Following  Seas  (0  degrees) 


Figure  A30.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  B, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  Following  Seas  (0  degrees) 
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Figure  A31.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  A, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  Head  Seas  (180  degrees) 


Figure  A32.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  B, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  Head  Seas  (180  degrees) 
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Figure  A33.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  A, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  Starboard  Beam  Seas  (90  degrees) 
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Figure  A34.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  B, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  Starboard  Beam  Seas  (90  degrees) 
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Figure  A3  5.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  A, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  Port  Beam  Seas  (270  degrees) 


Figure  A36.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  B, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  Port  Beam  Seas  (270  degrees) 
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Figure  A3  7.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  A, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  45  degrees  Seas. 


Figure  A38.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  B, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  45  degrees  Seas. 
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Figure  A3  9.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  A, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  135  degrees  Seas. 


Figure  A40.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  B, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  135  degrees  Seas. 
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Figure  A41.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  A, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  225  degrees  Seas. 


Figure  A42.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  B, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  225  degrees  Seas. 
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Figure  A43.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  A, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  315  degrees  Seas. 


Figure  A44.  Plot  of  Vertical  Motion  (^)  per  [ft]  Wave  Height  of  Connection  Point  B, 
Onboard  the  Ship  (CAPE-D)  and  the  Barge  (RRDF)  at  315  degrees  Seas. 
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APPENDIX  B 


Length  of  the  Ramp 

L  =  mft  (30.5m) 

Width  of  the  Ramp 

B  =  25 ft  (7.6m) 

Mass  Density  of  the  Ramp 

0.284  Ibfiin 

Young’s  Module  for  Stress  Calculations 

E  =  29x\0^psi  (200GPa) 

Spring  Stiffness  between  Ramp  and  Barge 

^  =  205.6x10'  Ibflft  (3xl0'[A/m]) 

Damping  between  Ramp  and  Barge 

c  =  4xl0'  Ibfslft  (5.9xlO'[A5/m]) 

Table  2  Parameters  used  in  the  Numerieal  Solution 


Figure  7.  Sketeh  of  the  Ship-Stem  Ramp-Barge  Arrangement. 
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Figure  B1 .  Effect  of  Ramp  on  Barge  (RRDF)  Heave  Amplitude  per  [ft]  Wave  Height, 

at  Following  Seas  (0  degrees). 


Figure  B2.  Effect  of  Ramp  on  Barge  (RRDF)  Heave  Amplitude  per  [ft]  Wave  Height, 

at  Quartering  Seas  (45  degrees). 
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Figure  B4. 
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Effect  of  Ramp  on  Barge  (RRDF)  Heave  Amplitude  per  [ft]  Wave  Height, 
at  Starboard  Beam  Seas  (90  degrees). 


Effect  of  Ramp  on  Barge  (RRDE)  Heave  Amplitude  per  [ft]  Wave  Height, 
at  Quartering  Seas  (135  degrees). 
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Figure  B5.  Effect  of  Ramp  on  Barge  (RRDF)  Heave  Amplitude  per  [ft]  Wave  Height, 

at  Head  Seas  (180  degrees). 


Figure  B6. 
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Effect  of  Ramp  on  Barge  (RRDF)  Heave  Amplitude  per  [ft]  Wave  Height, 
at  Quartering  Seas  (225  degrees). 
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Figure  B7.  Effect  of  Ramp  on  Barge  (RRDF)  Heave  Amplitude  per  [ft]  Wave  Height, 

at  Port  Beam  Seas  (270  degrees). 


Figure  B8.  Effect  of  Ramp  on  Barge  (RRDF)  Heave  Amplitude  per  [ft]  Wave  Height, 

at  Quartering  Seas  (315  degrees). 
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Figure  BIO. 


Effect  of  Ramp  on  Barge  (RRDF)  Roll  Amplitude  per  [ft]  Wave  Height,  at 
Following  Seas  (0  degrees). 


Effect  of  Ramp  on  Barge  (RRDF)  Roll  Amplitude  per  [ft]  Wave  Height,  at 
Quartering  Seas  (45  degrees). 
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Figure  B1 1.  Effect  of  Ramp  on  Barge  (RRDF)  Roll  Amplitude  per  [ft]  Wave  Height,  at 

Starboard  Beam  Seas  (90  degrees). 


Figure  B12.  Effect  of  Ramp  on  Barge  (RRDF)  Roll  Amplitude  per  [ft]  Wave  Height,  at 

Quartering  Seas  (135  degrees). 


83 


Figure  B13.  Effect  of  Ramp  on  Barge  (RRDF)  Roll  Amplitude  per  [ft]  Wave  Height,  at 

Head  Seas  (180  degrees). 


Figure  B14.  Effect  of  Ramp  on  Barge  (RRDF)  Roll  Amplitude  per  [ft]  Wave  Height,  at 

Quartering  Seas  (225  degrees). 
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Figure  B15.  Effect  of  Ramp  on  Barge  (RRDF)  Roll  Amplitude  per  [ft]  Wave  Height,  at 

Port  Beam  Seas  (270  degrees). 
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Figure  B16.  Effect  of  Ramp  on  Barge  (RRDF)  Roll  Amplitude  per  [ft]  Wave  Height,  at 

Quartering  Seas  (315  degrees). 
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Figure  B17.  Effect  of  Ramp  on  Barge  (RRDF)  Pitch  Amplitude  per  [ft]  Wave  Height, 

at  Following  Seas  (0  degrees). 


Figure  B18. 


Effect  of  Ramp  on  Barge  (RRDF)  Pitch  Amplitude  per  [ft]  Wave  Height, 
at  Quartering  Seas  (45  degrees). 
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Figure  B19.  Effect  of  Ramp  on  Barge  (RRDF)  Pitch  Amplitude  per  [ft]  Wave  Height, 

at  Starboard  Beam  Seas  (90  degrees). 


Figure  B20.  Effect  of  Ramp  on  Barge  (RRDE)  Pitch  Amplitude  per  [ft]  Wave  Height, 

at  Quartering  Seas  (135  degrees). 
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Figure  B21.  Effect  of  Ramp  on  Barge  (RRDF)  Pitch  Amplitude  per  [ft]  Wave  Height, 

at  Head  Seas  (180  degrees). 
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Figure  B22.  Effect  of  Ramp  on  Barge  (RRDF)  Pitch  Amplitude  per  [ft]  Wave  Height, 

at  Quartering  Seas  (225  degrees). 
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Figure  B23.  Effect  of  Ramp  on  Barge  (RRDF)  Pitch  Amplitude  per  [ft]  Wave  Height, 

at  Port  Beam  Seas  (270  degrees). 


Figure  B24.  Effect  of  Ramp  on  Barge  (RRDF)  Pitch  Amplitude  per  [ft]  Wave  Height, 

at  Quartering  Seas  (315  degrees). 
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Figure  B25.  Heave  Amplitude  per  [ft]  Wave  Height  for  the  Barge  (RRDF),  with  Ramp, 

at  Variable  Wave  Directions. 


Figure  B26.  Heave  Amplitude  per  [ft]  Wave  Height  for  the  Barge  (RRDF),  with  Ramp, 

at  Variable  Wave  Directions. 
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Figure  B27.  Roll  Amplitude  per  [ft]  Wave  Height  for  the  Barge  (RRDF),  with  Ramp, 

at  Variable  Wave  Direetions 


Figure  B28.  Roll  Amplitude  per  [ft]  Wave  Height  for  the  Barge  (RRDF),  with  Ramp, 

at  Variable  Wave  Direetions  . 
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Figure  B30. 


Pitch  Amplitude  per  [ft]  Wave  Height  for  the  Barge  (RRDF),  with  Ramp, 
at  Variable  Wave  Directions 


Pitch  Amplitude  per  [ft]  Wave  Height  for  the  Barge  (RRDF),  with  Ramp, 
at  Variable  Wave  Directions  . 
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Figure  B3 1 .  Amplitude  per  [ft]  Wave  Height  of  Relative  Motion  Between  Ramp’s  and 
Barge  Connection  Points  at  Following  Seas  (0  degrees) . 


Figure  B32.  Amplitude  per  [ft]  Wave  Height  of  Relative  Motion  Between  Ramp’s  and 
Barge  Connection  Points  at  Quartering  Seas  (45  degrees) . 
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Figure  B33.  Amplitude  per  [ft]  Wave  Height  of  Relative  Motion  Between  Ramp’s  and 
Barge  Connection  Points  at  Starboard  Beam  Seas  (90  degrees) . 
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Figure  B34.  Amplitude  per  [ft]  Wave  Height  of  Relative  Motion  Between  Ramp’s  and 
Barge  Connection  Points  at  Quartering  Seas  (135  degrees)  . 
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Figure  B35.  Amplitude  per  [ft]  Wave  Height  of  Relative  Motion  Between  Ramp’s  and 
Barge  Connection  Points  at  Head  Seas  (180  degrees) . 
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Figure  B36.  Amplitude  per  [ft]  Wave  Height  of  Relative  Motion  Between  Ramp’s  and 
Barge  Connection  Points  at  Quartering  Seas  (225  degrees)  . 
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Figure  B37.  Amplitude  per  [ft]  Wave  Height  of  Relative  Motion  Between  Ramp’s  and 
Barge  Connection  Points  at  Port  Beam  Seas  (270  degrees) . 


Figure  B38.  Amplitude  per  [ft]  Wave  Height  of  Relative  Motion  Between  Ramp’s  and 
Barge  Connection  Points  at  Quartering  Seas  (315  degrees),  . 
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Figure  B40. 


Period  [sec] 

Amplitude  per  [ft]  Wave  Height,  of  Relative  Motion  Between  Ramp’s  and 
Barge  Conneetion  Point  A,  for  Various  Headings.  . 


Period  [sec] 


Amplitude  per  [ft]  Wave  Height,  of  Relative  Motion  Between  Ramp’s  and 
Barge  Conneetion  Point  A,  for  Various  Headings.  . 
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B41 .  Twist  Angle  per  [ft]  Wave  Height  of  Ramp  at  Various  Headings  . 


Figure  B42.  Twist  Angle  per  [ft]  Wave  Height  of  Ramp  at  Various  Headings  . 
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45  Heading 


Figure  B43. 


Displacement  Amplitude  per  [ft]  Wave  Height  of  the  Ramp  at  Wave 
Period  18  sec  (0.35  Hz)  and  at  Various  Headings  . 


Figure  B44.  Displacement  Amplitude  per  [ft]  Wave  Height  of  the  Ramp  at  Wave 
Period  18  sec  (0.35  Hz)  and  at  Various  Headings  . 
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Figure  B45.  Displacement  Amplitude  per  [ft]  Wave  Height  of  the  Ramp  at  Wave 
Period  6.28  sec  (1  Hz)  and  at  Various  Headings  . 
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Figure  B46.  Displacement  Amplitude  per  [ft]  Wave  Height  of  the  Ramp  at  Wave 
Period  6.28  sec  (1  Hz)  and  at  Various  Headings  . 
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Figure  B47.  Ramp  Maximum  Stress  Amplitude  per  [ft]  Wave  Height  for  Various 

Headings  . 
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Figure  B48.  Ramp  Maximum  Stress  Amplitude  per  [ft]  Wave  Height  for  Various 

Headings  . 
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Figure  B49.  Ramp  Stress  Amplitude  per  [ft]  Wave  Fleight  at  Wave  Period  18  sec  (0.35 
Hz)  and  90  degrees  Heading  (k,  c  values  as  Figure  B43). 
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Figure  B5 1 .  Ramp  Stress  Amplitude  per  [ft]  Wave  Height  at  Wave  Period  6.28  sec  (1 
rad/sec)  and  90  degrees  Heading  (k,  c  values  as  Figure  B45). 


Figure  B52.  Ramp  Stress  Amplitude  per  [ft]  Wave  Height  at  Wave  Period  6.28  sec  (1 
rad/sec)  and  270  degrees  Heading  (k,  c  values  as  Figure  B46). 
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APPENDIX  C 


This  Appendix  contains  the  following  figures: 

•  Figure  Cl :  Pierson  -  Moskowitz  Spectrum. 

•  Figures  C2  through  C13:  Spectrum  of  the  relative  motion  response  at 
connection  point  A  for  all  cases  of  spring/damper  values  under 
consideration,  and  for  different  headings. 

•  Figures  C14  through  C20:  RMS  value  of  the  relative  motion  of  ramp  and 
barge  at  connection  point  A  for  all  cases  of  spring/damper  values 
considered.  Both  rectangular  and  polar  plots  are  presented  . 

•  Figures  C21  through  C27:  RMS  value  of  the  ramp  twist  for  all  cases  of 
spring/damper  values  considered  and  vertical  motions  at  the  Ramp/RRDF 
connection  points.  Both  rectangular  and  polar  plots  are  presented  here. 

•  Figures  C28  through  C34:  RMS  value  of  the  ramp  maximum  stress  for  all 
cases  of  spring/damper  values  considered  and  vertical  motions  at  the 
Ramp/RRDF  connection  points.  Both  rectangular  and  polar  plots  are 
presented  here. 

•  Figures  C35  through  C36:  RMS  value  of  the  vertical  motion  of  the  barge 
connection  point  A  and  B  for  all  cases  of  spring/damper  values 
considered. 

Discussion  of  these  figures  is  presented  in  Chapter  V,  Section  D. 
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Figure  Cl .  Pierson  -  Moskowitz  Speetrum 


Figure  C2.  Speetrum  of  the  Relative  Motion  Response  at  Conneetion  Point  A  at 
Following  Seas  (0  Degrees)  for  all  Cases  Under  Consideration. 
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Figure  C3.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  B  at 
Following  Seas  (0  degrees)  for  all  Cases  Under  Consideration. 


Figure  C4.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at 
Following  Seas  (0  degrees),  (holding  k  constant  and  varying  c). 
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Spectrum 


0.012 


Figure  C5.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at 
Following  Seas  (0  degrees),  (holding  k  constant  and  varying  c). 


Figure  C6.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at 
Following  Seas  (0  degrees),  (holding  k  constant  and  varying  c). 
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Spectrum 


Figure  C7.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at 
Quartering  Seas  (45  degrees)  for  all  Cases  Considered. 


Figure  C8.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at 
Starboard  Beam  Seas  (90  degrees)  for  all  Cases  Considered. 
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Figure  C9.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at 
Quartering  Seas  (135  degrees)  for  all  Cases  Considered. 


Figure  CIO.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at  Flead 
Seas  (180  degrees)  for  all  Cases  Considered. 
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Figure  C 1 1 .  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at 
Quartering  Seas  (225  degrees)  for  all  Cases  Considered. 


Figure  C12.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at  Port 
Beam  Seas  (270  degrees)  for  all  Cases  Considered. 
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Figure  C13.  Spectrum  of  the  Relative  Motion  Response  at  Connection  Point  A  at 
Quartering  Seas  (315  degrees)  for  all  Cases  Considered. 


Figure  C14.  RMS  Value  of  the  Relative  Motion  of  Ramp  and  Barge  at  Connection 

Point  A  for  all  Cases  Considered 
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Figure  C15. 


Figure  Cl 6. 


Polar  Plot  of  the  RMS  Value  of  the  Relative  Motion  of  Ramp  and  Barge  at 
Connection  Point  A  (holding  k  constant  and  varying  c). 


Polar  Plot  of  the  RMS  Value  of  the  Relative  Motion  of  Ramp  and  Barge  at 
Connection  Point  A  (holding  k  constant  and  varying  c). 
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Figure  Cl 7. 


Figure  C18. 


Polar  Plot  of  the  RMS  Value  of  the  Relative  Motion  of  Ramp  and  Barge  at 
Connection  Point  A  (holding  k  constant  and  varying  c). 


Polar  Plot  of  the  RMS  Value  of  the  Relative  Motion  of  Ramp  and  Barge  at 
Connection  Point  A  (holding  c  constant  and  varying  k). 
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Figure  Cl 9. 


Figure  C20. 


Polar  Plot  of  the  RMS  Value  of  the  Relative  Motion  of  Ramp  and  Barge  at 
Connection  Point  A  (holding  c  constant  and  varying  k). 


Polar  Plot  of  the  RMS  Value  of  the  Relative  Motion  of  Ramp  and  Barge  at 
Connection  Point  A  (holding  c  constant  and  varying  k). 
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RMS  Value  of  the  Ramp  Twist  for  all  Cases  Considered 


Figure  C21. 


Figure  C22.  Polar  Plot  of  the  RMS  Value  of  the  Ramp  Twist,  (holding  k  constant  and 

varying  c). 
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Figure  C23. 


Figure  C24. 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Twist,  (holding  k  eonstant  and 

varying  e). 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Twist,  (holding  k  constant  and 

varying  c). 
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Figure  C25. 


Figure  C26. 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Twist,  (holding  e  constant  and 

varying  k). 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Twist,  (holding  c  constant  and 

varying  k). 
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Figure  C27.  Polar  Plot  of  the  RMS  Value  of  the  Ramp  Twist,  (holding  c  constant  and 

varying  k). 


Figure  C28.  RMS  Value  of  the  Ramp  Maximum  Stress  for  all  the  Cases  Considered 
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Figure  C29. 


Figure  C30. 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Maximum  Stress,  (holding  k 
constant  and  varying  c). 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Maximum  Stress,  (holding  k 
constant  and  varying  c). 
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Figure  C31. 


Figure  C32. 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Maximum  Stress,  (holding  k 
constant  and  varying  c). 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Maximum  Stress,  (holding  c 
constant  and  varying  k). 
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Figure  C33. 


Figure  C34. 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Maximum  Stress,  (holding  c 
eonstant  and  varying  k). 


Polar  Plot  of  the  RMS  Value  of  the  Ramp  Maximum  Stress,  (holding  c 
eonstant  and  varying  k). 
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Figure  C35. 


Figure  C36. 
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Plot  of  the  RMS  Value  of  the  Barge  Vertieal  Motion  (^)  at  Conneetion 
Point  A  for  Cases  Considered. 
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Plot  of  the  RMS  Value  of  the  Barge  Vertieal  Motion  (^)  at  Conneetion 
Point  B  for  all  Cases  Considered. 
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